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ABSTRACT 


The  phenomenon  of  two-phase  slug  flow  in  the  system  air-water 
was  studied  in  a  one-half  inch  clear  Lucite  tube  40  feet  in  length. 
Pressure  drop,  pressure  fluctuations,  slug  velocity,  and  slug  frequency 
were  measured  over  the  18.29  foot  test  section.  These  variables  were 
studied  in  horizontal  and  inclined  flow,  at  plus  and  minus  seven  degrees 
to  the  horizontal.  The  data  have  been  cross-plotted  to  show  the  in¬ 
fluence  of  the  respective  mass  flow  rates  and  tube  orientation  on  an 
air-water  system  at  80  degrees  Fahrenheit  and  essentially  atmospheric 
pressure.  Comparisons  between  the  pressure  drop  correlations  of  Lockhart 
and  Marti nelli  and  of  Kordyban  were  made  with  the  data.  The  predictions 
were,  however,  found  to  be  between  30  to  50  percent  lower  than  the  actual 
data.  To  improve  on  this,  a  semi -empirical  theory  based  on  a  mechanistic 
flow  model  of  slug  flow  was  developed  and  compared  with  the  data.  The 
predictions  by  this  model  were  found  to  be  between  25  to  30  percent 
higher  than  the  data  over  most  of  the  range  studied.  By  providing  an 
upper  bound  to  slug  flow  pressure  drop,  this  model  is  superior  to 
either  of  the  above  two  models  for  slug  flow. 
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CHAPTER  I 


INTRODUCTION 

The  case  of  single  phase  flow  in  pipes  has  been  satisfactorily 
developed  and  it  is  conceded  that  the  flow  of  a  Newtonian  fluid  offers 
a  fairly  simple  problem.  However,  the  introduction  of  another  phase 
makes  the  task  considerably  more  complex.  It  might  seem  that  the  ad¬ 
dition  of  a  small  amount  of  liquid  to  a  gas  stream  would  still  allow 
the  pressure  drop  to  be  calculated  by  a  gas  flow  equation.  However, 
if  the  calculation  were  made  on  this  basis,  the  result  might  be  in 
error  by  as  much  as  five  hundred  percent.  Errors  of  this  order  of 
magnitude  indicate  that  there  is  little  understanding  of  the  basic 
mechanism  of  two-phase  flow. 

As  a  consequence  of  this  lack  of  knowledge  many  designers 
prefer  to  avoid  two-phase  flow  wherever  possible.  There  are  circum¬ 
stances,  however,  where  some  technical  or  economic  advantage  can  be 
gained  by  utilizing  the  characteristics  of  two-phase  flow,  such  as  in 
heat  transfer,  mass  transfer,  agitation,  or  transportation  of  multiple 
phases,  as  could  be  encountered  in  wet  gas  gathering  systems.  Con¬ 
sequently  some  predictability  of  the  fluid  behavior  must  be  available 
to  the  designer. 

In  general,  a  two-phase  mixture  separates  when  flowing  into 
a  pipe.  The  gas  flows  away  rapidly  and  leaves  the  liquid  to  move  along 
at  a  much  slower  velocity.  However,  when  the  fluids  encounter  ob¬ 
structions,  such  as  bends,  valves,  or  a  change  in  slope,  the  liquid  may 
bridge  the  tube  cross  section,  resulting  in  the  formation  of  slugs  of 
liquid,  separated  by  gas-filled  voids.  With  changes  in  fluid  properties. 


1 


?  ' 


2 


such  as  those  caused  by  the  expanding  gas,  and  changes  in  tube  charac¬ 
teristics  and  orientation,  the  flow  pattern  may  again  resolve  itself 
into  any  one  of  several  distinct  forms.  The  prediction  of  pressure  drop 
and  liquid  inventory  in  the  pipe  is  therefore  very  difficult.  For 
example,  Flanigan  (21)  reported  that  the  pressure  drop  in  a  pipeline 
had  been  predicted  to  be  30  psig.  This  soon  built  up  to  300  psig! 

It  is  therefore  evident  that  the  behavior  of  two-phase  flow 
cannot  be  predicted  by  the  utilization  of  the  existing  single  phase 
flow  equations.  The  fluid  mixture  cannot  be  simply  described  as  being 
laminar,  transitional,  or  turbulent.  In  single  phase  flow,  the  frictional 
pressure  drop  is  inversely  proportional  to  the  fifth  power  of  the  tube 
diameter.  In  multiphase  flow,  however,  this  diameter  can  be  effectively 
reduced  for  each  phase,  resulting  in  a  drastic  increase  in  pressure  drop. 
While  single  phase  momentum  transfer  occurs  by  either  a  molecular  or 
turbulent  mechanism,  two-phase  flow  has  an  additional  bulk  momentum 
transfer,  as  well  as  an  acceleration  term. 

It  is  remarkable  that  a  general,  but  empirically  derived, 
correlation  as  that  of  Lockhart  and  Martinelli  (30),  which  was  developed 
without  regard  for  flow  patterns,  can  be  applied  over  practically  the 
whole  range  of  variables.  But  neither  this,  nor  other  correlations 
developed  in  a  similar  manner,  or  from  a  separate  treatment  of  each 
phase,  adequately  describe  two-phase  flow.  If  one  knows  the  flow 
pattern  and  all  the  physical  and  geometric  properties  of  the  system, 
the  best  correlations  will  predict  the  pressure  drop  with  an  accuracy 
of  about  25  percent  (2,3).  To  overcome  this  problem  the  general  di¬ 
rection  of  research  has,  therefore,  been  directed,  in  recent  years, 
towards  obtaining  analytical  solutions  for  specific  flow  geometries. 
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The  most  frequently  cited  reference  for  flow  patterns  is 
Alves  (1).  He  categorized  the  patterns  into  seven  regions  based  on 
visual  observations  of  flow  in  horizontal  pipes.  In  order  of  their 
appearance  with  increasing  gas  rate  and  at  a  constant  liquid  rate 
they  are  illustrated  in  Figure  1  and  are  called: 

1 .  bubble  flow 

2.  plug  flow 

3.  stratified  flow 

4.  wavy  flow 

5.  slug  flow 

6.  annular  flow 

7.  dispersed  or  spray  flow. 

The  majority  of  these  flow  forms  have  been  studied  rather 
extensively,  with  the  exception  of  horizontal  slug  flow.  It  is  perhaps 
the  most  troublesome  form  of  flow  to  be  encountered  in  the  design  and 
operation  of  multiphase  pipelines. 

Slug  flow  is  formed  when  the  liquid  bridges  the  tube  cross 
section  and  is  therefore  characterized  by  a  pattern  in  which  all  or 
part  of  the  cross  section  of  the  tube  becomes  alternately  filled  with 
liquid  plugs,  separated  by  gas-filled  voids.  These  slugs  can  travel 
at  a  considerably  higher  velocity  than  would  be  expected  from  the 
overall  volumetric  rate.  When  vented  to  the  atmosphere,  slug  flow  is 
characteri zed  by  violent  pressure  fluctuations  which  are  propagated 
through  the  pipeline.  If  these  pulsations  are  not  damped,  they  could 
cause  severe  damage  to  the  receiving  vessels  at  the  end  of  the  line. 

To  reduce  the  effects  of  slug  flow,  or  to  completely  avoid  this 
troublesome  pattern  it  is  therefore  necessary  to  study  the  problem  in  detail. 
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FIG.  1  -  HORIZONTAL  FLOW  PATTERNS  AFTER  ALVES 
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It  was  not  until  1961  that  this  complex  flow  pattern  was  first  analyzed. 
Kordyban's  (27)  mechanistic  flow  model  for  slug  flow  was  used  to  corre¬ 
late  his  pressure  drop  data.  In  1963  he  (28)  experimentally  analyzed 
the  behavior  of  slug  flow  by  measuring  slug  velocities  and  frequencies, 
pressure  drops  and  fluctuations,  and  velocity  distributions.  The 
nature  of  this  work  was  largely  exploratory  since  it  was  carried  out 
in  the  slow  slug  flow  region.  The  work  was  also  left  incomplete. 

This  thesis  could,  perhaps,  be  considered  to  be  an  extension 
of  Kordyban's  work.  Similar  measurements  were  made,  but  in  the  region 
of  highly  slugging  flow.  In  addition,  slug  flow  in  pipe-lines  was 
simulated  and  studied  by  inclining  the  tube  to  positive  and  negative 
angles  of  seven  degrees  (20)  from  the  horizontal. 


CHAPTER  II 


LITERATURE  REVIEW 

Two-phase  flow  has  been  studied  periodically  for  the  past 
sixty  years,  but  it  was  not  until  the  early  1940 's  that  this  phenomenon 
was  begun  to  be  analyzed  in  detail. 

In  the  field  of  gas-liquid  flow  the  literature  is  quite 
extensive.  A  recent  annotated  bibliography  (26)  cites  over  2800  refer¬ 
ences  for  the  period  1950  to  1962.  Gouse  (23)  has  compiled  a  general 
index  of  5253  references  and  a  further  review  (12)  from  1931  to  1962 
lists  308. 

Initial  attempts  to  describe  the  behavior  of  a  gas-liquid 
mixture  in  terms  of  single  phase  flow  have  proved  to  be  largely  un¬ 
successful.  Frictional  pressure  drop  has  been  calculated  from  mixture 
fluid  properties  such  as  density  and  viscosity,  but  acceleration  effects 
due  to  the  expanding  gas  cannot  be  calculated  analogously.  The  arti¬ 
ficiality  of  the  assumption  of  homogeneous  flow,  and  particularly  of  the 
assumption  that  the  gas  and  liquid  velocities  are  the  same  has  been 
realized  for  some  time.  The  relative  movement,  or  slip,  between  the 
phases  has,  therefore,  been  studied  extensively.  This  phenomenon 
accounts  for  the  ratio  of  liquid  and  gas  in  the  pipe,  or  holdup,  which 
is  not  only  dependent  upon  the  relative  rates  of  the  phases  but  also  on 
tube  characteristics ,  orientation  and  fluid  properties. 

To  fully  understand  the  basic  mechanisms  of  two-phase  flow  it 
is  therefore  necessary  to  study  these  aspects  and  their  interrelation 
in  detail.  While  the  earliest  work  was  primarily  concerned  with  the 
study  and  prediction  of  pressure  drop,  increasing  emphasis  has  been 
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placed  on  analysis  of  the  individual  flow  patterns  encountered  in  gas- 
liquid  flow.  This  sequence,  however,  has  not  always  been  followed. 
Consequently  this  review  of  the  literature  has  been  divided  into  three 
sections:  pressure  drop,  flow  patterns,  and  holdup. 

2 . 1  Pressure  Drop 

The  earliest  two-phase  flow  literature  was  primarily  con¬ 
cerned  with  the  theory  and  application  of  air-lift  pumps.  A  discussion 
of  these  reports  by  Wood  (44)  points  out  that  a  simple  mathematical 
theory  on  the  action  of  the  airlift  pump  was  presented  by  Lorenz  (31) 
in  1909.  Davis  and  Weidner  (16)  in  1911  conducted  tests  with  vertical 
air-water  flow  and  determined  the  efficiency  of  actual  work  done  to 
that  obtainable  from  theory. 

In  1930  Uren  et  al  (42)  reported  on  the  behavior  of  air-oil 
mixtures.  He  derived  a  parameter  similar  to  a  Reynolds  number,  using 
a  specific  gravity  based  on  the  volume  fractions  of  the  air  and  oil, 
and  an  absolute  viscosity  calculated  from  Poiseuille's  formula.  A 
plot  of  Uren's  parameter  versus  the  Fanning  friction  factor,  on 
logarithmic  coordinates,  resulted  in  a  linear  relationship  for  all 
values  of  this  Reynolds  number.  It  was  realized,  however,  that  since 
slippage  and  non-homogeneity  of  the  air-oil  mixture  had  not  been  con¬ 
sidered  in  this  relationship,  it  was  open  to  criticism. 

A  study  by  Moore  and  Wilde  (35)  in  1931  took  slippage  into 
account.  A  comparison  of  the  actual  friction  drop  with  that  predicted 
by  the  Fanning  equation  by  subtracting  the  hydrostatic  and  velocity 
heads  from  the  total  pressure  drop  showed  large  deviations.  This 
procedure  was  therefore  considered  inapplicable  for  the  flow  of  hetero- 


geneous  fluids.  Nevertheless  their  empirical  expression  correlated 
their  data,  but  was  found  to  be  unsatisfactory  for  longer  test  sections 
such  as  in  gas-lift  oil  wells. 

The  analysis  of  data  from  28  producing  wells,  carried  out 
by  Nowels  (37)  in  1932,  resulted  in  a  "turbulence  factor"  analogous 
to  Reynolds  number.  The  parameters  were  calculated  using  an  average 
density  based  on  top  and  bottomhole  pressure  and  an  average  velocity 
calculated  from  the  volumetric  rates  of  the  oil  and  gas  under  average 
pressure.  The  absolute  viscosity  was  calculated  from  the  linear  re¬ 
lationship  of  oil-gas  mixture  viscosities  and  fraction  oil  in  the 
mixture  at  atmospheric  conditions.  Again,  as  in  the  work  of  Uren  et 
al  (42),  a  linear  relationship  was  obtained  in  the  logarithmic  plot 
of  "turbulence  factor"  versus  friction  factor.  The  result  was  that, 
for  a  given  value  of  "turbulence  factor",  the  friction  factor  increased 
with  the  diameter,  which  is  contradictory  to  the  effect  in  the  hori¬ 
zontal  flow  phenomenon.  Pressure  gradient  was  greatest  for  pipes  of 
small  diameter. 

Up  to  this  time  the  great  majority  of  publications  was  con¬ 
cerned  with  vertical  two-phase  flow.  In  1939,  Boelter  and  Kepner  (6) 
published  the  results  of  their  investigation  into  horizontal,  as  well 
as  inclined  flow.  Air-oil  and  air-water  mixtures  in  one-half  and  three- 
quarter  inch  standard  pipes  were  studied  in  the  horizontal  position 
and  at  an  upward  slope  of  1:6.  The  flow  was  observed  through  glass 
sections  and  a  variety  of  flow  patterns  were  described,  including  slug 
flow.  In  computing  the  "power"  [sic]  to  lift  oil  up  a  sloped  pipe  it 
was  assumed  that  the  oil  flowed  quiescently  in  the  lower  portion  of  the 
pipe.  The  "equilibrium  quantity"  [sic]  of  oil  within  the  pipe  was  de¬ 
termined  as  a  function  of  air-oil  ratio,  air  velocity,  and  pipe  diameter 
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The  authors  set  up  the  equipment  such  that  a  horizontal 
section  of  pipe  preceded  the  inclined  one.  Under  these  conditions 
the  pressure  gradient  in  the  sloped  line  was  observed  to  increase. 

In  addition  the  level  of  oil  on  this  section  increased,  thereby  causing 
the  oil  level  and  pressure  gradient  in  the  horizontal  section  to  increase 
as  well.  The  pressure  drop  in  the  sloped  line  was  found  to  be  almost 
doubled  at  high  oil  rates  and  10  cfm  air  rate,  while  the  increase  in 
the  horizontal  section  was  about  twenty-five  percent. 

In  1944,  Martinelli  and  co-workers  (34)  published  the  results 
of  a  study  of  isothermal  two-phase,  two-component  flow  of  air  and 
eight  liquids,  which  included  benzene,  water,  and  S.A.E.  40  oil.  The 
results  of  the  tests  which  were  conducted  in  a  one  inch  glass  pipe 
and  a  one-half  inch  galvanized  pipe  caused  the  authors  to  conclude 
that  tube  diameter,  liquid  and  gas  rates,  and  liquid  and  gas  viscosi¬ 
ties  and  densities  and  surface  tension  affected  pressure  drop.  Liquid 
viscosity  had  the  greatest  influence  on  pressure  drop  if  the  liquid 
was  flowing  "viscously"  [sic],  while  surface  tension  influenced  the 
flow  patterns  appreciably. 

In  developing  their  general  pressure  drop  correlation, 
Martinelli  et  al  (34)  assumed  that  the  two-phase  pressure  drop  is  a 
function  of  single  phase  pressure  drop  if  either  fluid  were  flowing 
alone.  Based  on  this  postulate,  these  four  possibilities  could  occur: 

1.  gas  turbulent  -  liquid  turbulent 

2.  gas  turbulent  -  liquid  viscous 

3.  gas  viscous  -  liquid  turbulent 

4.  gas  viscous  -  liquid  viscous 

The  first  three  mechanisms  were  analyzed,  although  only  the 


first  two  could  be  observed.  The  fourth,  it  was  stated,  probably 
could  not  occur-,  however,  Marti nelli  could  not  have  considered  very 
slow  stratified  flow,  where  laminar  flow  for  each  phase  is  a  possibility. 
In  the  development  of  their  correlations,  Martinelli  et  al  further 
postulated  that 

1.  The  flow  pattern  is  homogeneous  along  the  length  of  the 
tube.  This  implied  that  such  patterns  as  slug  and  plug  flow  were 
eliminated  from  consideration. 

2.  Pressure  drops  in  each  phase  are  equal. 

In  1949  Lockhart  and  Martinelli  (30),  using  additional  data 
with  tubes  from  0.0586  to  1.017  inches  inside  diameter  and  various 
liquids,  improved  on  the  correlations  of  Martinelli  et  al  (34).  A 
new  general  parameter  applicable  to  all  four  types  of  flow  was  pro¬ 
posed.  Furthermore,  slippage  between  the  phases,  determined  by  the 
fraction  of  the  tube  filled  with  each  phase,  i.e.,  holdup,  was  in¬ 
corporated. 

Alves  (1),  in  1954,  issued  his  report  on  the  study  of  two- 
phase  flow  for  the  design  of  pipe-line  contactors.  While  his  work 
was  primiarly  concerned  with  the  investigation  of  the  various  types 
of  flow  patterns  to  be  found  in  two-phase  flow,  he  also  reported 
pressure  drop  data.  His  slug  flow  pressure  drop  data  were  always 
higher  than  the  predictions  by  the  Lockhart-Martinel 1 i  correlation. 

In  1955  Chenoweth  and  Martin  (11)  proposed  an  improved 
correlation  for  high  pressure,  turbulent  two-phase  flow  in  horizontal 
pipes  of  large  diameter.  This  development  was  precipitated  by  the  in¬ 
applicability  of  the  Lockhart-Marti nel 1 i  (30)  correlation  to  predict 
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pressure  drop  for  pressures  above  50  psig  and  for  large  diameter  tubes. 

The  test  sizes  were  1.5  and  3.0  inches  tube  diameter.  The 
correlation  incorporated  a  parameter  which  was  equal  to  the  ratio  of 
the  fictitious  all-gas  to  all-liquid  system.  Furthermore,  any  fittings 
such  as  elbows,  valves,  and  orifices  as  may  be  found  in  the  test  section 
were  included. 

A  comparison  of  the  data  with  the  Lockhart-Martinel 1 i  corre¬ 
lation  showed  best  agreement  at  atmospheric  pressure.  The  largest  de¬ 
viations  were  found  to  occur  at  100  psia  in  the  3  inch  pipe,  where  the 
predicted  pressure  drop  was  higher  than  the  observed  by  a  factor  ranging 
from  1.4  to  2.5.  Chenoweth  and  Martin's  own  correlation,  on  the  other 
hand,  represented  all  of  the  data  within  plus  and  minus  50  percent, 
with  92  percent  of  the  data  falling  within  35  percent.  The  correlation 
was  also  found  to  fit  the  data  of  other  investigators,  including  those 
of  Alves  ( 1 ) . 

In  1957,  Hoogendoorn (24)  reported  his  results  on  air-water 
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and  air-oil  systems  in  horizontal  smooth  pipes  with  inner  diameters 
ranging  from  about  1.0  to  5.0  inches.  He  stated  that  the  Lockhart- 
Martinelli  (30)  correlation  was  not  valid  for  plug,  slug,  and  froth 
flow  if  the  gas  density  differed  from  that  of  air  at  atmospheric 
pressure.  For  wave,  mist,  and  annular  flow  it  was  inadequate  under 
any  conditions. 

Hoogendoorn  attempted  to  overcome  these  deficiences  by 
developing  a  new  correlation.  This  was  successful  only  for  the  lower 
range  of  gas  superficial  velocity.  The  standard  deviation  between  his 
data  and  the  correlation  was  6  percent.  In  the  higher  ranges,  where 
slug  flow  prevails,  the  standard  deviation  rose  to  14  percent.  By 
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comparison  with  Lockhart  Martinelli's  correlation,  the  standard  de¬ 
viation  from  his  data  was  13  percent,, 

The  same  year  Brigham  et  al  (7)  studied  gas-liquid  flow  in 
a  1.975  inch  pipe,  horizontal  and  inclined  at  12.4  degrees.  The  fluids 
used  were  air  and  No.  10  S.A.E.  oil  and  air  and  water.  The  data  were 
compared  with  the  correlation  of  White  (43).  For  the  regions  of  wave, 
cresting,  and  semi-annular  flow,  the  predicted  values  showed  good 
agreement  with  the  data.  In  slug  flow,  however,  the  measured  pressure 
drop  was  much  greater  than  predicted.  It  was  also  found  that  for 
nominally  stratified  flows  the  angle  of  slope  did  not  affect  the 
overall  pressure  drop.  In  a  comparison  of  the  inclined  slug  flow  data 
with  the  predictions  of  the  Lockhart-Martinel 1 i  (30)  equation,  the 
predicted  values  were  five  to  six  times  smaller  than  observed. 

With  the  data  he  collected  from  multiphase  pipeline  field 
tests,  Flanigan  (21)  in  1958,  modified  existing  design  equations  to 
improve  the  predictability  of  frictional  and  static  head  pressure 
drop.  His  analysis  showed  that 

1.  practically  all  of  the  pressure  drop  occurred  in  the  uphill 
sections  of  the  line  and, 

2.  the  pressure  drop  in  the  line  decreased  as  the  gas  rate  was 
i ncreased. 

The  latter  observation  is  contrary  to  that  observed  in  horizontal  flow, 
but  has  been  shown  to  be  true  in  vertical  lines  (10,  22,  39). 

The  pressure  drop  in  the  pipeline  was  obtained  by  applying 
a  pipeline  efficiency,  which  is  strongly  dependent  on  flow  patterns, 
to  the  frictional  pressure  drop.  The  static  head  was  also  corrected, 
since  the  two-phase  pipe  is  never  completely  filled  with  liquid.  Com- 


pari  son  of  predicted  with  actual  pressure  drop  was  found  to  result  in 
a  maximum  deviation  of  ±  14  percent,  with  the  majority  falling  within 
±  5  percent. 

One  of  the  first  mechanistic  flow  models  for  slug  flow  was 
presented  by  Kordyban  (27)  in  1961.  The  resulting  equation  resembles 
rather  closely  the  applicable  portion  of  the  empirical  Chenoweth- 
Martin  (11)  correlation.  The  two-phase  pressure  drop  was  calculated, 
in  each  case,  by  dividing  the  single  phase  pressure  drop  by  the  liquid 
volume  fraction  raised  to  a  power  between  0.75  and  1.00.  Kordyban' s 
correlation  was  developed  by  considering  all  the  liquid  to  flow  in 
slugs,  separated  by  gas-filled  spaces,  with  the  gas  pressure  drop 
assumed  to  be  negligible.  Only  frictional  pressure  drop  was  considered. 
However,  he  realized  that  in  the  steam-water  system  the  pressure  drop 
produced  by  acceleration  of  liquid  due  to  increasing  amounts  of  vapor 
could  be  calculated  to  be  over  50  percent  of  the  total  pressure  drop. 

For  the  0.420  inch  tube  agreement  of  the  prediction  with  his  data  was 
quite  good;  however,  for  the  0.315  inch  tube  the  calculated  values 
were  found  to  be  20  to  30  percent  lower  than  the  experimental  data. 

In  1963  Kordyban  and  Ranov  (28)  issued  a  further  report  on 
slug  flow  in  horizontal  tubes.  Pressure  drop  in  the  1.25  inch  plexi- 
glas  tube  was  measured  with  differential  pressure  transducers  rather 
than  with  U-tube  manometers,  as  had  been  the  usual  practice  up  to 
this  time.  This  enabled  them  to  measure  not  only  average  values,  but 
also  the  pressure  fluctuations  which  occur  when  slugging  flow  is 
vented  to  the  atmosphere.  At  the  time  of  its  presentation  the  work 
was  incomplete,  thus,  no  new  correlation  was  proposed.  The  approach 
in  this  study  of  slug  flow  was  to  measure  the  characteristic  behavior 
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of  the  slug  mechanism.  Variables  such  as  slug  velocity  and  frequency, 
pressure  drop  and  fluctuations,  and  velocity  distribution  were  there¬ 
fore  measured  to  enable  the  authors  to  improve  on  the  previous  corre¬ 
lation  of  Kordyban  (27) .  This  earlier  mathematical  model,  the  authors 
stated,  "did  not  describe  satisfactorily  the  behavior  of  slug  flow, 
because  it  had  been  based  on  insufficient  data". 

In  1964  a  critical  review  of  the  then  known  correlations  for 
two-phase  flow  was  reported  by  Dukler  and  Wicks  (17).  A  search  of 
the  literature  revealed  over  20,000  individual  measurements  of  pressure 
drop.  By  statistically  culling  the  available  data,  2620  points  re¬ 
mained.  An  overall  comparison  of  five  correlations  with  the  data 
showed  the  Lockhart-Marti nel 1 i  (30)  correlation  to  be  clearly  superior 
to  those  of  Baker  (4),  Bankoff  (5),  Chenoweth-Marti n  (11),  and  Yagi  (45) 

Dukler  et  al  (18)  further  investigated  frictional  pressure 
drop  through  similarity  analysis.  The  basic  equation  was  solved  for 
four  sets  of  boundary  conditions.  Comparison  of  these  results  to  the 
Lockhart-Marti nel 1 i  (30)  correlation  showed  them,  in  general,  to  be 
superior  to  the  latter,  while  the  opposite  was  true  at  low  viscosities 
with  small  diameter  tubes. 

Hughmark  (25)  in  1965  presented  a  correlation  which  utilized 
the  amount  of  liquid  in  situ,  i „ e . ,  holdup,  in  horizontal  slug  flow. 

This  holdup  was  estimated  from  the  relationship  of  bubble  velocity  and 
liquid  slug  Reynolds  number.  Comparison  with  the  results  of  several 
other  investigators,  including  Alves  (1),  Baker  (4),  and  Brigham  (7) 
was  inconclusive.  The  authors  stated  that  this  occurred  "probably  be¬ 
cause  some  of  the  latter  data  were  not  for  slug  flow." 

Chisholm  (13),  in  1966,  presented  a  note  on  the  relationships 
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between  friction  and  liquid  cross  section  during  two-phase  horizontal 
flow.  He  correlated  data  to  within  20  percent  in  smooth  and  galvanized 
tubes . 

2.2  Flow  Patterns 

The  first  reported  flow  pattern  observations  were  those  of 
Davis  and  Weidner  (16)  in  1911.  They  found  that  with  an  increasing 
air  rate  at  a  constant  water  rate  the  flow  pattern  in  a  vertical  tube 
would  change  from  one  of  small  bubbles  in  the  water  to  long  bubbles 
or  pistons,  which  disappeared  when  the  air  to  water  ratio  became  large. 

Boelter  and  Kepner  (6),  who  in  1939  carried  out  studies 
of  both  horizontal  and  inclined  flow,  installed  glass  sections  into 
their  1/2  inch  and  3/4  inch  tubes  for  flow  pattern  observations. 

Various  patterns  were  noted,  among  them  slug  flow.  The  difference  in 
behavior  of  oil  and  water  under  the  influence  of  a  high  air  rate  was 
also  described.  In  the  transition  from  wavy  flow,  the  wave  crests  of 
oil  flow  tended  to  break  up  into  spray,  while  water  formed  a  slug  flow 
pattern. 

A  year  later,  in  1940,  Cromer  and  Huntington  (15)  conducted 
extensive  investigations  using  2-inch  standard  pipe  fitted  with  ob¬ 
servation  windows.  Patterns  of  bubble,  slug,  froth,  annular  and  dis¬ 
persed  flow  were  observed  with  increasing  flow  rate  of  air. 

In  1944,  Martinelli  et  al  (34),  gave  a  detailed  description 
of  the  change  in  flow  patterns  with  a  decrease  in  water  rate  and  a 
constant  gas  rate.  Their  conclusions  were  that  "the  microscopic  be¬ 
havior  of  the  two-phase  flow  system  is  extremely  complex.  Since, 
however,  the  static  pressure  drops  from  the  water  and  water  with  surface 
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tension  reducing  additive  were  almost  identical,  although  the  detailed 
flow  patterns  were  quite  different,  it  appears  that  it  should  be  possi¬ 
ble  to  analyze  the  problem  from  the  macroscopic  point  of  view  and  ob¬ 
tain  results  which  are  of  utility  from  an  engineering  standpoint,  even 
though  the  details  of  flow  mechanism  are  not  clarified."  This  conclusion, 
the  authors  felt,  was  justified  in  the  same  manner  as  the  employment 
of  a  friction  factor  for  the  prediction  of  pressure  drop  in  turbulent 
flow,  without  the  need  for  consideration  of  the  complex  details  of 
turbulence.  As  a  result  a  two-phase  flow  system  was  considered  to 
contain  the  two  fluids  as  a  homogeneous  mixture. 

Lockhart  and  Martinelli  (30)  reconfirmed  these  conclusions 
in  1949.  Only  homogeneous  flow,  i.e„,  complete  mixing  of  the  phases, 
was  considered  in  their  correlations . 

In  the  discussion  of  this  paper  (30)  Carl  Gazeley,  Jr.  and 
O.P.  Bergelin  (University  of  Delaware  at  Newark,  Delaware)  concluded, 
however,  that  a  completely  successful  correlation  for  two-phase  flow 
could  not  be  obtained  until  the  transition  points  between  the  various 
types  of  flow  had  been  determined. 

The  same  year  Kosterin  (29)  published  his  results  on  the 
influence  of  tube  diameter  and  position  upon  hydraulic  resistance  and 
flow  structure  of  gas-liquid  mixtures.  Visual  observations,  still 
photographs,  and  motion  pictures  of  gas  liquid  flow  structures  in 
tubes  of  1,  2,  3,  and  4  inches  were  used  to  construct  charts  showing 
the  distribution  of  the  various  regions.  Slug  flow  was  described  as 
occurring  in  the  region  between  divided  flow  and  strongly  dispersed 
flow,  where  the  gas  or  vapor  moved  in  the  form  of  large  bubbles,  which 
occupied  a  considerable  part  of  the  diameter. 


The  possible  flow  configurations  were  divided  into  five 
distinct  regions. 

lc  Divided  flow:  with  gas  in  the  upper  part  and  liquid  in 
the  lower  part  of  the  tube  cross  section. 

2.  Quiet  plug  flow:  with  large  gas-bubble  plugs,  with  a  clean 
boundary  surface  and  without  froth  formation. 

3.  Plug  flow:  with  froth  formation  in  the  trailing  part  of 
the  bubble. 

4.  Plug  flow:  with  froth  formation  over  the  entire  gas-liquid 
boundary  surface  along  the  gas  plugs. 

5.  Strongly  dispersed  flow,  in  which  the  gas-liquid  mixture 
flows  in  the  form  of  a  more  or  less  uniform  froth:  a  gas-liquid 
emulsion. 

Flow  condition  charts  were  constructed  with  percent  gas  as 
the  ordinate  and  mixture  velocity,  based  on  the  gravimetric  delivery 
rate  of  the  mixture,  as  the  abscissa. 

In  inclined  flow,  at  70  degrees  to  the  horizontal,  it  was 
noted  that  stratified  flow  did  not  occur  and  thus,  it  was  stated,  was 
a  phenomenon  particular  to  horizontal  flow,  only. 

In  1954,  Alves  (1),  in  carrying  out  his  work  on  cocurrent 
pipeline  contactors,  gave  an  even  more  detailed  description  of  the 
various  possible  flow  configurations.  He  classified  horizontal  two- 
phase  flow  into  seven  more  or  less  distinctive  regions  (Figure  1). 

With  an  increasing  gas  rate  and  at  a  constant  liquid  rate,  they  are  (8) 

1.  Bubble  flow:  Bubbles  of  gas  move  along  the  upper  part  of 
the  pipe  at  approximately  the  same  velocity  as  the  liquid. 

2.  Plug  flow:  Alternate  plugs  of  liquid  and  gas  move  along  the 
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upper  part  of  the  pipe* 

3*  Stratified  flow:  Liquid  flows  along  the  bottom  of  the  pipe 
and  gas  flows  above,  over  a  smooth  liquid-gas  interface. 

4.  Wavy  flow:  Similar  to  stratified  flow  except  that  the  gas 
moves  at  a  higher  velocity  and  the  interface  is  disturbed  by  waves 
travelling  in  the  direction  of  flow. 

5.  Slug  flow:  Waves  are  picked  up  by  the  more  rapidly  moving 
gas  to  form  a  slug,  which  passes  through  the  pipe  at  a  velocity  greater 
than  the  average  liquid  rate. 

6.  Annular  flow:  The  liquid  flows  in  a  thin  film  around  the 
inside  wall  of  the  pipe  and  the  gas  flows  at  a  high  velocity  as  a 
central  core. 

7.  Dispersed  or  spray  flow:  Most  of  the  liquid  is  entrained 
as  a  spray  by  the  gas. 

A  plot  depicting  these  regions  of  two-phase  flow,  shown  in 
Figure  2,  was  set  up  by  Baker  (4)  in  the  same  year.  Valid  flow  pattern 
predictions  can  frequently  be  made  by  this  plot,  but  more  often,  only 
an  indication  of  the  possible  flow  regime  can  be  expected  for  a  given 
set  of  conditions. 

In  1957,  Hoogendoorn  (24),  constructed  flow  pattern  charts 
using  the  parameters  of  Kosterin  (29).  The  results  were  similar,  but 
because  the  only  variables  taken  into  account  were  mass  velocity  and 
percent  gas,  the  charts,  in  each  case,  were  particular  only  to  the 
specific  system  under  investigation. 

Several  such  area  plots  have  been  constructed  in  the  course 
of  investigations  into  two-phase  flow.  However,  researchers  in  the 
field  generally  agree  (8)  that,  to  date,  there  is  no  method  in  the 
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FIG.  2  -  BAKER'S  CHART  FOR  TWO-PHASE  FLOW  REGIONS 
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published  literature  which  adequately  describes  or  establishes  just 
what  type  of  flow  will  occur  under  the  range  of  possible  conditions. 

2.3  Holdup  and  Slippage 

The  apparent  slippage  of  the  gas  through  the  liquid  was 
realized  by  the  earliest  investigators,  including  Lorenz  (31)  in  the 
study  of  vertical  two-phase  flow  systems.  The  first  reported  measure¬ 
ments  of  holdup  were  those  of  Moore  and  Wilde  (35)  in  1931.  Quick 
closing  valves  were  employed  to  trap  the  mixture.  The  authors  re¬ 
ported  that  holdup  was  independent  of  viscosity. 

For  horizontal  flow,  the  first  reported  measurements  and 
correlation  of  holdup  were  those  of  Lockhart  and  Martinelli  (30). 

The  liquid  was  trapped  in  the  tube,  then  washed  out  with  large  amounts 
of  volatile  solvent,  and  the  excess  solvent  was  evaporated.  A  plot 
of  holdup  versus  the  square  root  of  the  ratio  of  liquid  to  gas  pressure 
drop  resulted,  for  all  practical  purposes,  in  a  single  curve. 

In  1957  Hoogendoorn  (24)  measured  holdup  by  means  of  a 
capacitive  probe.  His  results  were  checked  by  adding  an  oil-soluble 
radioactive  tracer  to  measure  the  radiation  intensity,  which  depends 
linearly  on  liquid  holdup.  In  a  comparison  of  Lockhart-Martinel 1 i 1 s 
(30)  correlation  with  these  data,  Hoogendoorn  reported  deviations ,  but 
their  severity  was  not  disclosed.  He  correlated  his  own  data  with  the 
type  of  slip  velocity  relationship  normally  employed  in  vertical  flow 
(35).  No  significant  effect  of  pipe  diameter  or  liquid  viscosity  was 
found  in  these  velocity  ranges. 

The  adverse  effects  of  holdup  in  pipelines  are  clearly  evident 
from  the  report  of  Flanigan  (21)  in  1957.  As  had  been  reported  pre- 
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viously,  the  predicted  pressure  drop  of  30  psig  in  a  pipeline  soon  rose 
to  300  psig.  This  drastic  increase  was  attributed  to  holdup  occurring 
in  the  line. 

In  1963  an  improved  method  for  holdup  determination  was  used 
by  Neal  and  Bankoff  (36),  who  employed  an  electrical  resistivity  probe. 
The  circuit  consisted  of  a  steel  needle  connected  to  a  flip-flop  circuit. 
The  needle  was  traversed  across  the  pipe  diameter  to  measure  bubble 
frequency  at  various  positions.  This  procedure  was  found  to  be  a  great 
improvement  on  the  earlier  method  in  which  the  probe  was  stationary. 
Petrick  (38)  also  used  the  traversing  technique  but  found  that  it  still 
only  gave  an  average  value. 

Kordyban  (28)  in  1963,  used  a  resistance  probe  to  measure 
holdup.  The  apparatus  utilized  the  variation  of  electrical  resistance 
with  liquid  level  between  two  parallel  electrodes  which  extended  almost 
completely  across  the  tube  diameter.  Although  slugs  and  waves  could 
be  detected,  these  results  were  affected  by  a  lag  in  drainage  of  the 
liquid  film  from  the  probe. 

In  1963  McManus  (33)  reviewed  experimental  techniques  in 
two-phase  flow.  In  the  section  on  the  measurement  on  density  or  void 
fractions,  reference  was  made  to  the  development  of  contact  probes,  con¬ 
ductance  elements,  radiation  attenuation,  and  light  absorption.  His 
conclusion  was  that  these  methods  will  only  give  overall  density  or 
void  fraction  at  pipe  cross  sections.  Precise  determination  of  mean 
film  depth  and  surface  character  in  stratified  and  annular  flows  are 
possible,  however,  non-homogeneous  flows  cannot  be  measured  successfully. 

Other  investigators  such  as  Hughmark  (25)  continued  to  employ 
methods  of  estimating  holdup  by  means  of  ratio  of  the  respective  phase 
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volumetric  rate  to  total  volumetric  rate,  corrected  with  an  experi¬ 
mentally  determined  factor. 

The  main  purpose  of  the  above  methods  of  determining  holdup 
was  to  develop  correlations  for  predicting  holdup  and  consequently 
pressure  drop.  Most  of  these  methods  are  inadequate,  as  no  completely 
reliable  instrument  for  measuring  holdup  is  available. 

As  a  practical  solution  to  the  problem  of  holdup  in  multi¬ 
phase  pipelines,  McDonald  and  Baker  (32)  developed  an  improved  method 
of  "pigging".  Whereas  pigging  was  tradi tional ly  used  to  remove 
hydrates,  etc.,  this  new  method  employs  inflatable  rubber  spheres  to 
remove  the  accumulated  liquid.  The  use  of  these  pigs,  introduced  into 
the  line  at  predetermi ned  intervals,  was  found  to  drastically  increase 
pipeline  capacity.  Although  this  method  is  in  wide  scale  use,  it  is 
economical  only  if  the  ratio  of  liquid  to  gas  is  relatively  small.  As 
the  relative  amount  of  liquid  is  increased,  the  authors  stated,  the 
advantages  from  sphere  operation  will  decrease.  Furthermore,  the 
optimum  interval  between  the  introduction  of  pigs  cannot  be  determined 
precisely  and  thus  the  rate  of  pigging  is  based  largely  on  the  pre¬ 
determined  maximum  allowable  pressure  drop  in  the  line. 


CHAPTER  III 


THEORY  OF  SLUG  FLOW 

Two  of  the  existing  two-phase  flow  correlations  have  been 
employed  in  this  thesis  as  a  comparison  with  the  data.  The  corre¬ 
lation  of  Lockhart-Marti nel 1 i  (30)  was  empirically  derived  for  a  homo¬ 
geneous  flow  pattern,  whereas  that  of  Kordyban  (26)  was  developed 
from  a  simplified  model  of  slug  flow.  Although  the  first  is  widely 
used,  a  pattern  such  as  slug  flow  was  ruled  out  by  the  assumptions 
made  in  the  derivation.  Kordyban's  model  was  found  to  be  oversimplified 
since  it  was  "based  on  insufficient  data"  (28). 

Visual  observations  and  the  measurements  of  slug  flow  carried 
out  by  Kordyban  in  his  later  report  (28)  suggested  the  model  shown 
below.  It  was  observed  that  the  liquid  film  at  the  front  of  a  slug 
is  picked  up  by  the  slug  and  its  velocity  is  increased.  A  similar 
amount  of  liquid  is  dropped  off  at  the  rear  of  this  slug  and  its  mo¬ 
mentum  is  gradually  lost  until  the  film  is  again  picked  up  by  the 
following  slug. 


3.1  Slug  Flow  Model 
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Consider  a  slug  moving  at  a  constant  velocity  us„  The  film, 
which  is  flowing  at  velocity  u,  is  picked  up  by  the  slug  and  is  dropped 
off  at  the  rear.  It  is  assumed  that  the  film  has  been  accelerated  to 
the  velocity  of  the  slug,  u  .  If  it  is  further  assumed  that  the  slug 
does  not  accelerate,  a  momentum  balance  may  be  written  around  the  control 
surface.  This  results  in  the  following  equation: 


ZF  =  0  =  Pus^S  -  p<u2>  S  +  2ttR 
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Since  the  film  is  picked  up  by  each  slug  in  the  tube,  the  total  pressure 
drop  is  proportional  to  the  number  of  slugs  in  the  tube,  . 
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The  shear  force  on  the  slug  can  be  approximated  by  the  Fanning 
Equation  ;  therefore 
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Kordyban  (27)  regarded  only  this  shear  face  term  to  be  sig¬ 
nificant.  He  calculated  the  velocity  of  the  slugs  by  assuming  the 
liquid  to  flow  alone  in  the  tube. 

To  improve  on  this,  the  shear  face  term  was  evaluated  in  this 
thesis  by  assuming  the  gas  to  be  incompressible.  Thus  the  slug  velocity 
equals  the  total  fluid  velocity, 


Since  the  length  of  the  slugs  is  not  known,  the  fraction  of 
the  tube  occupied  by  liquid  slugs  can  be  approximated  as 

NsLs  ■  L  <^> 


where  L$  =  length  of  slugs 

L  =  length  of  test  section. 

Therefore  equation  (4)  becomes 
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To  evaluate  the  momentum  exchange  term  in  equation  (3) 


N  pS 


AP  =  -4—  (u  -  <u2>) 
m  A  v  s  ' 


(5) 


(6) 


several  simplifying  assumptions  must  be  made.  The  film  cross  sectional 
area,  S,  may  be  calculated  from  the  holdup,  R^,  which  was  measured  by 
Lockhart  and  Martinelli  (30)  by  trapping  the  liquid  in  the  tube.  The 
liquid  in  the  film  is  therefore  equal  to  the  total  liquid  less  the 
liquid  due  to  the  flow  rate. 
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Therefore 


S  =  A  ( R,  - 


l  QL  +  QG 


)  =  Ac 


and  equation  (6)  becomes, 

AP  =  N  pe(u  2  -  <u2>)  (7) 

mss 

2  2 

For  simplicity,  <u  >  -  <u>  .  The  resulting  error  of  16.5%,  if  the 
film  is  assumed  to  be  in  laminar  flow,  is  justified  in  view  of  more 
inexact  assumptions  made  above. 

Evaluation  of  <u>  is  difficult  since  the  film  flows  in  the 
bottom  of  the  tube.  However,  if  it  is  assumed  that  the  film  flows  on 
a  flat  plate,  an  approximation  of  <u>  can  be  found  by  employing  the 
Navier  Stokes  Equation. 
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The  equation  of  continuity  is 


3u  9u  3u 
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To  solve  these  equations,  several  simplifying  assumptions  must 
be  made.  Since  the  flow  is  only  in  the  x  direction,  u^  and  uz  are  zero. 
Therefore  the  continuity  equation  reduces  to 
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and  therefore 
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and  there  is  no  change  in  u  with  respect  to  z 
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The  gas  phase  pressure  gradient  on  the  film  in  the  x  di¬ 
rection  is  negligible,  therefore 


The  equation  is  therefore  reduced  to 


which  must  be  solved  for  u  . 

A 

This  form,  however,  is  similar  to  the  heat  transfer  equati 
for  temperature  in  a  slab. 
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A  solution  to  the  film  slowdown  problem  may  therefore  be 
found  by  considering  the  momentum  loss  to  occur  analogously  to  heat 
loss . 


Again  considering  the  diagram, 
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the  velocity  of  the  film  leaving  the  slug  is  equal  to  the  slug  velocity 
u$.  The  film  then  begins  to  slow  down  since  the  velocity  at  the  wall 
is  equal  to  zero.  To  determine  the  average  film  velocity  in  the  x 
direction,  <u>,  it  is  necessary  to  solve  the  partial  differential 
equati on. 

A  solution  to  the  equation,  with  respect  to  y  and  t  is  (14) 
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For  the  film,  u(y,0)  =  Uq  =  ug. 

To  find  the  average  velocity  in  the  film, 
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k  =  y/p 

t  =  time  interval  between  slugs 
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-  inverse  of  slug  frequency 
l  =  depth  of  liquid  in  the  film 
r  fraction  liquid  in  film  to  tube  diameter 


=  (R,  - 


L  \  + 


)  D  =  eD 


Substitution  into  the  approximated  change  in  momentum 
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yields  the  momentum  loss  term  of  the  total  pressure  drop 


o£  ,  82u  2  -  .  .{il ndiVktj 

APm  =  Ns  fc  '  -7s-  — —*  e  H  >  > 


(2n-l  )' 


or 


-  N,  ££  u.2  {1  -  M 

TT 


S  g.  s 
3c 


00 


({ 


(2n-l  )' 


r(2n-l)Vktn 

e  4£2  )2} 


To  this  term  is  added  the  friction  loss  of  the  slugs 

2fpu  A  Q, 
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where  f  is  approximated  by 
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The  momentum  term  is  divided  by  g  to  maintain  consistent 
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The  overall  pressure  drop  equation  for  the  slug  flow  model 


therefore  becomes 
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A  sample  calculation  is  shown  in  Appendix  C. 


CHAPTER  IV 


EXPERIMENTAL  EQUIPMENT 

The  experimental  apparatus,  as  shown  in  the  schematic  flow 
diagram.  Figure  3,  and  the  photographs,  Figure  4,  consisted  of  two 
main  parts: 

(a)  metering  system 

(b)  pipeline. 

The  design  and  construction  of  the  equipment  was  based 
largely  on  findings  obtained  from  preliminary  investigations  carried 
out  on  bench  scale  equipment.  This  apparatus  consisted  of  a  0.59 
inch  inside  diameter  glass  tube,  about  8  feet  long,  into  which  water 
and  air  from  the  laboratory  bench  were  introduced  via  a  mixing  tee. 
Sizing  was  carried  out  on  the  basis  of  Baker's  flow  regime  plot  (4). 
Aside  from  using  the  initial  test  set-up  to  study  fluid  behavior,  it 
was  also  used  for  the  development  of  an  electronic  slug  sensor  de¬ 
signed  to  trace  out  the  profile  of  the  individual  slugs  for  possible 
holdup  correlation  (40). 

The  test  section,  consisting  of  clear  Lucite  tubing,  of 
one-half  inch  diameter,  had  a  total  length  of  40  feet.  The  seven 
sections  were  joined  by  flanges,  sealed  with  rubber  0-rings.  Since 
no  data  on  entrance  length  in  slug  flow  were  reported  in  the  liter¬ 
ature,  the  first  pressure  tap  was  located  almost  halfway  along  the 
tube,  17.9  feet  downstream  from  the  mixing  tee.  This  gave  an  en¬ 
trance  length  of  430  pipe  diameters.  The  second  pressure  tap  was 
located  18.29  feet  further  downstream,  resulting  in  an  exit  length 
of  4.0  feet.  The  pipeline  was  mounted  on  the  edge  of  three  laminations 


31 


32 


FIG.  3  SCHEMATIC  FLOW  DIAGRAM  OF  THE  EQUIPMENT 
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METERING  SYSTEM 


PIPELINE 


FIG.  4  -  EXPERIMENTAL  EQUIPMENT 
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of  3/4  inch  plywood,  and  the  whole  apparatus  was  suspended  on  a  ful¬ 
crum,  which  was  located  about  halfway  along  the  length  for  ease  in 
altering  the  inclination  of  the  tube.  For  added  rigidity,  the  whole 
assembly  was  clamped  to  rectangular  frames  fastened  to  the  wall. 

In  this  apparatus  the  water  was  circulated  in  a  closed  loop, 
while  the  air  was  obtained  from  the  compressed  air  system  and  was 
vented  to  the  atmosphere  after  passing  through  the  test  section. 

The  water  circulation  was  carried  out  from  a  2  foot  diameter,  45 
gallon,  sheet  iron  tank  equipped  with  a  3/8  inch  copper  tube  cooling 
coil  and  two  100  watt  Chromalox  heaters (Model  RI-100).  The  water  was 
pumped  from  this  tank  by  a  1/2  horse  power  Jabsco  pump  (Model  JC-1/2"- 
37)  with  1/2  inch  ports.  The  outlet  stream  was  split  up  with  part 
going  through  a  bypass,  part  being  used  to  saturate  the  preheated 
incoming  airstream,  and  the  remainder  going  to  the  flow  meters.  The 
flowrate  was  measured  by  two  Brooks  rotameters  (Model  No.  8-1100  and 
Model  No.  6-111 0-GTM ) .  The  metered  water  stream  was  conducted  to  the 
mixing  tee,  at  the  far  end  of  the  test  section,  through  3/4  inch 
copper  tubing.  The  temperature  was  controlled  at  80  ±1  degrees  Fahren¬ 
heit  in  the  collection  tank  by  means  of  the  cooling  coil  and  heaters. 

The  cooling  water  was  circulated  at  a  constant  rate,  while  heat  input 
to  the  system  was  controlled  with  a  Chromalox  thermostat  (Model  AR  1014) 
and  a  Powerstat  variable  autotransformer  (Type  116B). 

To  saturate  the  entering  air  at  80  degrees  Fahrenheit,  it  was 
passed  through  a  steam  heated  Heat-X  (Model  120)  heat  exchanger  where 
its  temperature  was  raised  to  178  degrees  Fahrenheit  for  adiabatic  satu¬ 
ration  at  80  degrees  (9).  Quenching  to  80  degrees  was  carried  out  in 
a  special  tank  into  which  part  of  the  circulating  water  stream  was 
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sprayed  through  a  nozzle.  The  saturated  air  then  passed  through  a 
liquid  trap  and  to  two  Brooks  Rotameters (Model  No.  8-1100  and  Model 
No.  6-1110-GTM).  The  metered  air  was  then  passed  to  the  mixing  tee 
through  1/2  inch  copper  tubing.  The  flow  to  the  various  flow  meters 
and  tanks,  downstream  of  the  pump,  was  controlled  with  Crane  1/2  inch 
needle  valves  (Cat.  88).  The  two  fluids  were  introduced  to  the  test 
section  through  a  mixing  tee,  a  chamber  about  1-1/4  inches  long  and 
3/4  inches  in  diameter,  with  the  air  and  water  entering  at  right 
angles  to  each  other.  Two  125  mesh  stainless  steel  screens  were 
placed  downstream  from  the  mixing  tee  for  improved  mixing  and  to 
cause  a  pressure  drop  to  reduce  the  effect  of  pressure  surges  in  the 
test  section.  Pressure  surges  and  backflow  of  the  water  were  also 
controlled  by  the  surge  tank  located  ahead  of  the  mixing  tee. 

Pressure  drop  and  slug  presence  were  determined  by  two 
Statham  (Model  No.  PM80TC  ±15-350,  rated  for  ±15  psid  at  5000  psig) 
differential  pressure  transducers,  connected  to  the  test  section  by 
means  of  short,  water  filled,  3/16  inch  Polyflo  tubing.  Excitation 
of  these  transducers  and  amplification  of  the  signal  from  them  was 
carried  out  by  Hewlett-Packard  (Sanborn  Model  8805A)  carrier  pre¬ 
amplifiers,  operated  according  to  the  instruction  manual  (8805A).  The 
two  signals  were  recorded  simultaneously  on  a  Hewlett-Packard  (Moseley 
Autograf  Model  7100  B)  two  channel  stripchart  recorder. 

Slug  frequency  and  velocity  were  measured  by  various  means. 
These  data  were  primarily  taken  from  the  strip  charts,  however,  visual 
and  electronic  checking  was  done  over  most  of  the  range.  Low  slug 
frequencies  and  velocities  were  measured  visually  with  the  aid  of  a 
stopwatch.  At  higher  frequencies  the  slugs  were  sensed  with  a  photo- 
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conductivity  cell  (Clairex  CL5M7)  and  counted  on  a  Beckman  (Model  7350R) 
Universal  EPUT  and  Timer. 

After  passing  through  the  pipeline,  the  air  and  water  were 
allowed  to  separate  with  the  air  venting  to  the  atmosphere,  and  the 
water  returning  to  the  circulation  tank.  To  reduce  the  effect  of 
the  pressure  surges,  caused  by  the  slugging,  the  separator  tank  was 
packed  with  wire  mesh  screen  to  dissipate  the  energy  of  the  slugs. 


CHAPTER  V 


EXPERIMENTAL  PROCEDURE 

After  a  steady  temperature  of  80  degrees  Fahrenheit  was 
attained  in  the  circulation  tank,  the  air  was  turned  on,  as  well  as 
the  steam  to  the  air  heater.  The  heater's  temperature  was  adjusted 
to  180  degrees  by  manipulation  of  the  inlet  and  outlet  throttling 
valve.  Then  the  pump  was  switched  on,  allowing  the  water  to  regain 
steady  temperature  by  circulation. 

In  the  meantime  the  electronic  apparatus  was  allowed  to 
warm  up,  calibration  of  the  preamplifiers  was  checked,  and  the  strip 
chart  recorder  was  zeroed.  When  steady  temperatures  were  regained 
throughout  the  system,  the  valves  to  the  air  rotameters  were  opened 
first.  Then  the  water  rotameters  were  opened.  This  order  prevented 
flooding  of  the  air  system  with  water.  Finally,  the  water  line  to 
the  spray  tank  was  opened  and  the  needle  valve  at  the  bottom  of  this 
tank  was  coarsely  adjusted  to  maintain  a  liquid  trap  in  this  vessel. 
Testing  was  ready  to  begin. 

Data  were  taken  by  holding  the  air  rate  at  a  constant  level, 
and  advancing  the  water  rate  through  the  region  of  slug  flow.  This 
procedure  had  several  advantages  over  holding  the  water  rate  steadily 
and  advancing  the  air  rate.  It  was  found,  as  may  be  confirmed  from 
Baker's  chart,  Figure  2,  or  from  the  appearance  of  flow  patterns  ac¬ 
cording  to  Alves  (1)  that  if  the  water  rate  were  held  constant,  and 
the  air  rate  varied,  several  flow  patterns  had  to  be  passed  through 
before  the  slug  flow  pattern  was  reached.  With  a  steady  air  rate  and 
varying  water  rate,  however,  only  stratified  or  wavy  flow  had  to  be 
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passed.  An  additional  reason  for  using  the  preferred  procedure  was 
that  the  needle  valves  to  the  rotameters  were  not  precision  made  and 
were  therefore  affected  by  vibration  and  pressure  on  the  stems,  causing 
the  valves  to  close  gradually  and  therefore  required  almost  constant 
valve  manipulation  to  maintain  a  steady  water  rate.  This  effect  was 
not  noted  with  the  valves  to  the  air  rotameters,  although  the  same 
type  of  valve  was  employed. 

The  pressure  levels  and  fluctuations  from  each  of  the  two 
pressure  transducers  were  recorded  on  the  two-channel  strip  chart 
recorder.  Checks  of  slug  velocity  and  frequency  were  made  by  timing 
the  passage  of  a  slug  between  the  pressure  taps  and  by  counting  the 
number  of  slugs  in  a  certain  time  interval,  respectively.  Slug  fre¬ 
quency  was  further  determined  by  measuring  the  change  in  resistance 
in  a  photoconductivity  cell,  mounted  on  the  test  tubing,  when  a  slug 
of  liquid  passed  in  front  of  it.  To  amplify  the  change  in  resistance 
the  water  in  the  system  was  dyed  with  Flourescein  dye.  The  signal 
was  filtered  for  line  frequency,  and  the  slugs,  represented  by  oscil¬ 
lations,  were  counted  on  the  Beckman  EPUT  Meter. 
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CHAPTER  VI 


INTERPRETATION  OF  RESULTS 


6.1  Flow  Patterns 

The  transition  between  flow  patterns  at  a  constant  liquid 
rate  with  increasing  gas  flow  rate  was  observed  to  occur  as  described 
by  Alves  (1)  and  as  is  shown  on  Baker's  chart,  Figure  2.  With  an  in¬ 
creasing  air  rate  and  a  low  water  rate,  the  order  of  patterns  observed 
was  stratified,  wavy  and  annular.  At  a  somewhat  higher  water  rate, 
stratified,  plug,  slug  and  annular  flow  occurred.  At  the  high  liquid 
rates  the  order  of  bubble,  slug,  and  annular  flow  was  followed.  Mist 
flow  was  not  observed  due  to  pressure  limitations  of  the  compressed  air 
system. 

No  clear  distinction  appears  to  exist  between  plug  and  slug 
flow,  as  evidenced  by  the  dotted  transition  line  on  Baker's  area  plot, 
Figure  2.  At  low  fluid  rates  the  water  tended  to  bridge  the  tube  cross 
section,  forming  a  fairly  regular  series  of  long,  quietly  moving  water 
plugs.  With  an  increased  air  rate  the  rate  of  slug  formation  became 
somewhat  more  erratic  and  the  slug  length  became  quite  variable.  Froth 
began  to  form  at  the  front  of  the  slugs.  At  the  highest  air  rates  the 
whole  slug  consisted  of  a  frothy  mass  but  was  still  recognizable  as  a 
slug  by  its  stability  and  high  velocity.  A  further  increase  in  air 
rate  caused  the  collapse  of  slugs  with  the  subsequent  beginnings  of 
semi-annular  flow. 

The  range  of  the  slug  flow  pattern  that  was  studied  in  this 
report  has  been  superimposed  on  the  flow  pattern  map  of  Baker  (4). 

The  regions  covered  by  each  of  the  three  tube  inclinations  are  shown 
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in  Figures  5  to  7.  The  region  in  which  irregular  slug  flow  occurred, 
i.e.,  slugs  were  formed  infrequently  and  often  broke  up  soon  after 
formation,  was  marked  as  "Intermittent  Slugging"  on  these  plots.  This 
area  was  also  found  to  be  located  largely  outside  of  the  slug  flow 
region  denoted  by  Baker  (4).  It  is  suspected  that  the  entrance  condi¬ 
tions  are  responsible  for  this  phenomenon.  Hoogendoorn  (24)  noted  that 
by  using  different  entrance  conditions  it  was  possible  to  obtain  either 
slug  or  wave  flow  the  other  conditions  remaining  the  same.  Furthermore, 
Brigham  (7)  stated  that  he  was  unable  to  produce  stratified  flow  due 
to  the  return  bends  which  caused  the  liquid  to  bridge  the  tube  cross 
section.  A  further  observation  made  by  this  author  was  that  in  es¬ 
sentially  wavy  flow  the  slugs  leaving  the  outlet  of  the  pipeline  ap¬ 
peared  to  generate  new  slugs  at  the  mixing  tee.  It  was  also  possible 
to  create  slug  flow  by  momentarily  plugging  the  exit  to  cause  a  pressure 
buildup.  The  sudden  release  in  pressure  with  removal  of  the  obstruction 
caused  the  liquid  to  bridge  and  slugs  were  generated.  Scott  (41)  sum¬ 
marized  that  there  is  a  possibility  of  two  distinctly  different  flow 
patterns  occurring,  depending  on  such  factors  as  the  inlet  arrangement. 
Although  one  of  these  patterns  may  be  unstable,  its  transition  to  the 
more  stable  configuration  may  not  be  observable  within  the  limited 
scale  of  an  experimental  apparatus.  McDonald  (32)  also  verified  this 
point  in  describing  the  effect  of  water  and  air  being  fed  into  a  75 
foot  long  experimental  pipeline  of  1.25  inches  diameter.  The  flow 
pattern  was  observed  to  change  from  stratified  to  wave,  and  then  slug 
flow.  While  this  is  thought  to  be  primarily  the  effect  of  the  increasing 
gas  volumetric  fraction  due  to  the  gas  expansion  along  the  pipeline  it 
shows  that  the  flow  pattern  can  be  highly  unstable. 
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FIG.  5  -  COMPARISON  OF  BAKER'S  FLOW  PATTERN  CORRELATION  WITH  THE  DATA  IN  HORIZONTAL  FLOW 
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FIG.  6  -  COMPARISON  OF  BAKER'S  FLOW  PATTERN  CORRELATION  WITH  THE  DATA  In  POSITIVELY  SLOPED  FLOW  (+7 
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FIG.  7  -  COMPARISON  OF  BAKER'S  FLOW  PATTERN  CORRELATION  WITH  THE  DATA  IN  NEGATIVELY  SLOPED  FLOW 
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In  inclined  flow  Kosterin  (29)  noted  that  the  entrance  condi¬ 
tions  exerted  only  a  very  small  influence  on  the  areas  of  stable  flow 
structures.  Figure  6,  for  positively  inclined  flow,  shows,  however, 
that  the  slug  flow  area  is  even  further  enlarged  than  was  the  case  for 
horizontal  flow,  Figure  5.  The  reason  for  these  contradictory  obser¬ 
vations  is  probably  due  to  the  fact  that  Kosterin' s  tube  was  inclined 
at  70  degrees  to  the  horizontal,  virtually  vertical  flow,  while  the 
pipe  of  this  work  was  inclined  to  only  7  degrees.  The  intermittent 
slugging  noted  in  horizontal  flow,  which  was  primarily  caused  by  slug 
breakup,  was  still  evident.  No  stratified  flow  patterns  could  be  ob¬ 
served  in  this  tube  inclination,  as  confirmed  by  Brigham  (7)  and 
Kosterin  (29).  It  is  reasonable  to  conclude  that  gravity  effects 
preclude  stratified  flow  in  positively  inclined  tubes. 

In  downward  sloping  flow  the  region  of  intermittent  slugging 
was  found  to  be  almost  non-existent,  primarily  caused  by  slug  breakup 
occurring  before  the  test  section  was  reached.  The  region  of  strati¬ 
fied  and  wavy  flow  became  greatly  enlarged  and  slug  flow  did  not  occur 
until  the  higher  mass  rates  were  employed.  Under  the  circumstances  of 
negatively  inclined  tubes,  it  can  be  concluded  that  the  separated  flow 
patterns  would  be  favored  at  low  fluid  rates.  It  must  be  emphasized, 
however,  that  this  observation  may  not  be  true  for  tubes  inclined  at 
angles  greater  than  minus  seven  degrees. 

A  comparison  of  the  transition  lines  between  intermittent  and 
regular  slug  flow  for  all  three  inclinations  shows  them  to  be  virtually 
coincident.  To  forestall  the  argument  that  these  lines  may  have  been 
subject  to  operator  bias,  measurement  of  the  range  of  slug  flow  was 
carried  out  by  two  separate  means,  and  identical  areas,  for  both  types 


of  slug  flow,  were  obtained.  In  the  first  investigation  the  water  rate 
was  held  constant  while  the  air  rate  was  varied  for  each  setting  of 
water  rate.  In  the  second  run  the  air  was  held  constant  and  the  water 
rate  was  advanced,  again  through  the  complete  range  of  both  fluids. 

Each  test  run  was  carried  out  independently  of  the  other  and  the  plots 
of  the  two  sets  of  data  were  found  to  coincide.  It  must  therefore  be 
concluded  that  in  regularly  slugging  flow,  the  angle  of  inclination 
of  plus  and  minus  seven  degrees  does  not  alter  appreciably  the  area 
for  stable  slug  flow,  if  shown  on  Baker's  (4)  plot.  This  observation 
has  not  been  reported  in  the  literature  for  slug  flow.  However, 

Brigham  (7)  found  that  "the  angle  of  slope  definitely  does  not  affect 
the  overall  pressure  drop  when  the  fluids  are  flowing  in  wave,  cresting 
and  semi-annular  flow".  This  would,  of  course,  neglect  the  effect  of 
static  head.  Since  pressure  drop  is  directly  related  to  flow  patterns, 
regularly  slugging  flow  must  be  included  in  this  grouping  of  Brigham's. 

The  actual  data  for  this  work  were  taken  only  in  the  regular 
slug  flow  region  as  the  area  of  intermittent  slugging  may  have  been 
the  result  of  the  equipment  characteristics, as  discussed  above. 

6.2  Slug  Flow  Characteristics 

A  typical  example  of  a  strip  chart  recording  is  shown  in 
Figure  8.  The  pertinent  reference  lines  have  been  added  to  show  how 
the  measurement  for  the  data  points  was  carried  out.  Pressure  drop, 
pressure  fluctuations,  slug  frequency  and  slug  velocity,  which  was 
calculated  from  the  average  time  taken  for  a  slug  to  travel  the  length 
of  the  test  section,  were  read  off  the  charts.  Checks  on  velocity  and 
frequency  were  made  visually  and  electronically. 
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FIG.  8  -  STRIPCHART  RECORDING 
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6.2.1  Overall  Pressure  Drop 

The  measurements  of  total  pressure  drop  are  shown  in  Figures 
9  to  11 .  In  comparing  these  curves  with  respect  to  tube  inclination 
it  must  be  noted  that  the  inclined  tubes  were  not  corrected  for  static 
head,  which  is  a  function  of  holdup  (21).  As  a  result,  only  the  re¬ 
lative  magnitudes  and  slopes  of  the  individual  curves  are  discussed 
below. 

The  figures  show  the  strong  influence  of  gas  rate  on  overall 
pressure  drop.  At  the  lowest  air  rates  the  curves  are  found  to  be  quite 
flat  and  an  increasing  liquid  rate  causes  only  a  gradual  increase  in 
pressure  drop.  An  increased  air  rate,  however,  causes  the  slope  of  the 
curves  to  grow  much  more  rapidly.  There  is  an  upper  limit,  however. 

As  the  maximum  rates  were  attained,  the  slug  flow  pattern  changed  to 
semi-annular  and  annular  flow  since  slug  breakup  began  to  occur.  The 
effect  of  this  phenomenon  is  shown  in  the  gradual  decrease  in  slope  of 
the  upper  curves.  It  is  unfortunate  that  equipment  limitations  pre¬ 
vented  further  data  points  from  being  taken. 

A  comparison  of  the  results  of  the  three  tube  inclinations 
shows  the  curves  to  be  very  similar  with  respect  to  location  on  the 
chart  and  to  slopes.  The  offsets  are  due  primarily  to  the  effect  of 
static  head,  however,  they  may  also  be  influenced  by  the  lag  in  move¬ 
ment  of  the  recorder  pens  as  well  as  by  operator  bias.  The  difficulty 
in  taking  measurements  off  the  charts  lies  with  the  often  extremely 
violent  pressure  fluctuations  characteristic  of  slug  flow.  The  recorder 
pen  lag  would  have  reduced  this  effect  somewhat,  and  therefore  could 
have  caused  the  readings  to  be  in  error  by  several  percent.  This 
problem  could  be  overcome  in  part  through  the  use  of  a  recorder  with 
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FIG.  9  -  OVERALL  PRESSURE  DROP  FOR  HORIZONTAL  FLOW 
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..10  -  OVERALL  PRESSURE  DROP  FOR  POSITIVELY  SLOPED  FLOW  (+7°) 
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FIG.  11  -  OVERALL  PRESSURE  DROP  FOR  NEGATIVELY  SLOPED  FLOW  (-7°) 
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a  faster  response,, 

6.2.2  Pressure  Fluctuations 

Figures  12  to  14  show  the  magnitude  of  slug  flow  pressure 
fluctuations.  Although  this  phenomenon  has  been  noted  by  almost  all 
investigators  with  reference  to  slug  flow,  Kordyban  (28)  was  the  first 
to  publish  data  on  such  measurements.  The  use  of  differential  pressure 
transducers  makes  it  possible  to  record  the  fluctuations,  while  this 
was  not  possible  with  the  previously  employed  U-tube  manometers. 

Brigham  (7)  noted  that  it  was  necessary  to  design  special  buffers  in 
the  manometer  lines  to  reduce  this  effect. 

The  figures  show  that  at  low  gas  rates  with  increasing  water 
rate  the  pulsations  are  virtually  constant.  Again,  as  was  the  case 
for  overall  pressure  drop,  the  gas  rate  shows  a  greater  influence  than 
does  the  liquid  rate.  The  slopes  of  the  curves  increase  sharply  with 
increased  gas  rates,  but  again  the  gradual  change  to  semi-annular  flow 
is  noticeable.  With  a  constant  air  rate  and  increasing  water  rate  a 
maximum  is  obtained,  after  which  the  fluctuations  are  reduced  and  tend 
to  level  off.  The  pressure  surges  of  the  maximum  air  rates  are  actually 
found  to  be  lower  than  those  of  the  next  to  highest  air  rates.  The 
transition  from  very  slow  slug  flow,  commonly  referred  to  as  plug  flow, 
through  normal  slug  flow  and  to  semi-annular  flow  can  be  seen  on  these 
figures  by  the  changes  from  very  low  fluctuations  through  the  maximum 
and  the  final  decrease. 

No  plausible  explanation  could  be  found  for  the  fact  that 
the  maxima  in  horizontal  flow  were  reached  at  much  higher  water  rates 
than  for  either  of  the  two  tube  inclinations. 

The  gross  magnitude  of  these  pulsations  has  never  been  re- 
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FIG.  12  -  PRESSURE  FLUCTUATIONS  FOR  HORIZONTAL  FLOW 
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FIG.  13  -  PRESSURE  FLUCTUATIONS  FOR  POSITIVELY  SLOPED  FLOW  (+7°) 
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ported  in  the  literature,  although  it  is  normal  practice  for  multi¬ 
phase  gathering  systems  to  employ  a  buffer  at  the  terminals  to  prevent 
damage  to  the  receiving  vessels.  These  excessive  pressure  surges 
could,  however,  also  be  found  within  the  pipeline.  In  this  investigation 
it  was  necessary  to  pack  the  separator  tank  with  wire  mesh  screens  to 
dissipate  the  energy  of  the  slugs, which  caused  severe  vibrations  to  the 
equipment  at  high  rates.  Figures  12  to  14  show  that  these  surges  were 
in  the  order  of  25  to  30  percent  of  the  overall  pressure  drop. 

A  tentative  explanation  of  the  relative  magnitude  of  these 
pressure  fluctuation  is  presented  in  Figure  15.  Four  views  of  the 
motion  of  slugs  within  a  pipe  are  shown  with  respect  to  time.  In  the 
initial  situation  three  slugs  are  in  the  pipe.  The  second  sequence 
shows  the  slugs  to  have  advanced  to  the  right  by  one-half  slug  length 
and  the  third  and  fourth  view  shows  them  to  have  advanced  a  further 
two-half  lengths.  A  new  slug  appears  at  the  left  in  sequence  4.  In 
the  second  view  one-half  of  the  right-most  slug  has  already  left  the 
tube.  To  consider  the  pressure  profile  along  the  tube,  two  main  simpli¬ 
fying  assumptions  were  made: 

1.  Acceleration  effects  are  negligible 

2.  Pressure  drop  due  to  the  gas  is  negligible 

Each  of  the  slugs  in  the  tube  causes,  in  this  case  only,  a 
constant  pressure  drop  across  each  slug.  In  the  lower  part  of  Figure 
15,  then,  path  number  1  is  followed  from  point  A.  Each  slug  represents 
a  certain  pressure  drop,  while  the  gas-filled  voids  are  represented 
as  a  constant  pressure  line  from  the  front  of  the  left  slug  to  the  tail 
of  the  next.  The  gauge  pressure  goes  to  zero  at  the  outlet.  In  the 
second  view  of  the  pipe  half  of  the  right  slug  has  already  left  the 
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tube.  The  pressure  B  at  the  left  is  therefore  lower  by  an  amount  equal 
to  the  pressure  drop  of  the  half  slug  that  has  left  and  path  2  is  followed. 
The  sequence  is  continued  through  view  number  3  in  which  the  pressure 
at  the  left  is  at  C,  lower  than  A  by  one  whole  slug's  pressure  drop. 

In  view  4  the  pressure  has  again  risen  to  point  B  since  a  new  half  slug 
has  appeared  at  the  left.  This  whole  sequence  from  1  to  3  repeats  it¬ 
self  then  as  long  as  new  slugs  appear  at  the  left.  Therefore  the 
pressure  fluctuations  cycle  from  a  maximum  of  A  through  the  minimum  of 
C. 

The  magnitude  of  these  pressure  surges,  if  the  above  assumptions 
are  still  valid,  would  be  approximately  constant,  whether  the  outlet 
pressure  is  at  zero  or  at  a  positive  pressure.  The  differential  pressure 
would  therefore  be  constant  for  either  case,  while,  the  pressure  surges, 
relative  to  the  outlet  pressure,  could  be  negligible  with  respect  to  a 
high  exit  pressure.  For  a  low  outlet  pressure,  as  in  this  work,  they 
could,  however,  be  quite  substantial,  i.e.,  in  the  order  of  25  to  30 
percent  of  the  pressure  drop. 

No  attempt  has  been  made  to  correlate  the  data  with  this 
model.  It  also  remains  to  be  shown  what  effect  different  fluid  pro¬ 
perties  and  tube  sizes  would  have  on  the  relative  magnitude  of  these 
fluctuations. 

6.2.3  Slug  Velocity 

The  slug  velocity  curves,  shown  in  Figures  16  to  18,  indicate 
that  the  angle  of  tube  inclination,  within  ±7  degrees  from  the  hori¬ 
zontal,  has  little  effect  on  either  the  shape  of  the  curves  or  on  the 
magnitude  of  the  velocities.  As  was  the  case  for  pressure  drop  at  low 
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FIG.  16  -  SLUG  VELOCITY  FOR  HORIZONTAL  FLOW 
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FIG.  17  -  SLUG  VELOCITY  FOR  POSITIVELY  SLOPED  FLOW  (+7°) 
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air  rates,  an  increasing  water  rate  has  a  lesser  effect  on  the  change  in 
magnitude  then  does  the  air  rate.  It  is  found  that  the  slopes  of  the 
curves  rise  sharply  with  an  increasing  air  rate.  Again,  as  in  the  case 
for  pressure  drop,  a  maximum  slope  is  reached  with  the  highest  fluid 
rates  after  which  a  decrease  occurs.  This  effect  is  also  attributed 
to  the  transition  from  slug  to  semi-annular  flow. 

The  velocities  shown  on  the  graphs  represent  an  average  of 
the  velocities  of  individual  slugs.  Velocity  was  found  to  be  quite  vari¬ 
able  and  appeared  to  be  strongly  influenced  by  slug  size,  as  would  be 
expected  by  a  material  balance.  Although  both  methods  for  measuring 
the  average  time  required  for  a  slug  to  travel  the  length  of  the  test 
section-direct  timing,  and  reading  off  the  strip  charts  -  agreed  fairly 
well,  a  distinct  scatter  of  the  data  points  can  be  seen  on  the  velocity 
curves.  Nevertheless,  the  general  shape  of  the  curves  is  similar  for 
all  three  tube  inclinations.  Comparison  with  Kordyban's  (28)  results 
is  difficult  because  his  work  was  carried  out  at  lower  fluid  velocities, 
but  the  trend  is  again  very  similar.  The  air  rate  also  appeared  to 
have  a  greater  influence  on  slug  velocity  than  did  the  water  rate. 

Kordyban  (27,  28)  was  not  able  to  determine  acceleration  of 
the  slugs  since  his  work  was  carried  out  at  essentially  atmospheric 
pressure.  The  pressure  drop  in  this  work,  however,  rose  to  a  maximum 
of  over  six  psig,  so  that  acceleration  effects  would  certainly  not  be 
negligible.  By  taking  the  effect  of  the  liquid  film  in  reducing  the 
available  cross  sectional  area  for  flow  into  consideration,  the  measured 
slug  velocity  was  still  found  to  be  up  to  50  percent  higher  than  the 
calculated  fluid  velocity.  But  slug  velocity  is  directly  proportional 
to  the  pressure  differential  across  the  slug  so  that  the  gas  expansion, 
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caused  by  an  overall  pressure  drop  of  over  6  p s i g s could  accelerate  a 
slug  by  over  25  percent.  Slug  acceleration  was,  however,  not  measured 
in  this  work. 

A  further  effect  on  increased  velocity  could  be  the  effect 
of  the  froth  at  higher  rates.  The  froth  was  observed  to  increase  the 
thickness  of  the  film  as  well  as  to  increase  the  size  of  the  slugs  so 
that,  at  higher  rates,  the  film  may  no  longer  be  completely  picked  up 
by  the  slugs.  This  would  cause  a  restriction  in  the  tube  diameter 
and  result  in  higher  slug  velocities.  The  possibility  that  the  hold¬ 
up  predictions  by  Lockhart-Martinel 1 i 1 s  (30)  correlation  are  too  low 
could  also  be  a  cause  for  the  apparent  discrepancy  between  calculated 
and  measured  velocities.  It  would  appear,  however,  that  slug  acceler¬ 
ation  would  cause  the  major  contribution.  Kordyban  (27)  noted  that 
acceleration  due  to  increasing  amounts  of  vapor  in  his  steam-water 
system  could  be  calculated  to  be  over  50  percent  of  the  total  pressure 
drop,  even  though  his  work  was  carried  out  at  essentially  atmospheric 
pressure . 

6.2.4  Slug  Frequency 

The  effect  of  fluid  rates  and  tube  inclination  on  slug  fre¬ 
quency  is  shown  in  Figures  19  to  21.  At  low  gas  and  liquid  rates  the 
frequency  is  low.  It  was  observed  that  the  slugs  were  very  long,  about 
12  to  18  inches,  and  moved  at  a  low  velocity.  As  the  air  rate  was  in¬ 
creased,  the  frequency  did  not  change  appreciably.  At  very  low  air 
rates  and  an  increasing  water  rate,  however,  the  frequency  was  greatly 
increased,  while  a  steadily  increasing  air  rate  caused  the  frequency 
to  drop  and  then  level  off  for  each  water  rate. 
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MASS  RATE  OF  AIR  (LB/HR) 


FIG.  20  -  SLUG  FREQUENCY  FOR  POSITIVELY  SLOPED  FLOW  (+7°) 
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MASS  RATE  OF  AIR  (LB/HR) 

FIG.  21  -  SLUG  FREQUENCY  FOR  NEGATIVELY  SLOPED  FLOW  (-7°) 
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The  change  in  slope  in  the  constant  water  rate  curves  is  be¬ 
lieved  to  be  due  to  a  flow  pattern  transition.  It  appears  that  the 
change  in  slope  represents  the  transition  from  plug  flow  to  slug  flow 
since  it  coincides  very  closely  with  the  change  from  the  flow  pattern 
of  clean,  froth-free,  plugs  to  one  in  which  there  is  an  ever  increasing 
amount  of  froth  within  the  slug.  The  possibility  that  this  phenomenon 
is  a  function  of  entrance  conditions  should,  however,  not  be  ruled  out. 
Kordyban's  (28)  work,  although  carried  out  in  a  much  lower  region  of 
slug  flow,  does  show  a  similar  occurrence  for  his  data.  At  his  maximum 
water  rate  the  shape  of  the  curve,  plotted  on  the  same  axes  as  in  this 
work,  is  almost  identical  with  those  of  this  work.  It  could  therefore 
be  concluded  that  the  effect  of  entrance  conditions  is  negligible  with 
respect  to  the  shape  of  the  frequency  curves,  even  though  Kordyban's 
entrance  section  consisted  of  a  ten  inch  section  of  vertical  tubing, 
while  that  of  this  work  had  17.9  ft.  of  straight  tubing. 

The  only  effect  of  tube  inclination  on  slug  frequency  appears 
to  be  the  location  of  the  discontinuity.  In  horizontal  flow  the  change 
is  much  more  gradual  than  for  the  other  two  slopes,  with  not  much  dif¬ 
ference  being  apparent  between  the  latter  two. 

6 . 3  Comparison  of  the  Pressure  Drop  Correlations  with  the  Data 

Figures  22  to  27  show  a  comparison  of  the  pressure  drop  corre¬ 
lations  of  Lockhart-Martinel 1 i  (30)  and  Kordyban  (27),  and  the  slug 
theory  developed  in  this  work,  with  the  data.  Although  the  data  from 
all  three  tube  inclinations  have  been  presented,  static  head  was  not 
considered  in  the  correlations.  As  a  result,  comparisons  can  only  be 
made  with  the  horizontal  data. 
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CORRELATIONS 

1  SLUG 

2  LOCKHART-MARTINELLI 

3  KORDYBAN 
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a  POSITIVE  (+7°) 


FIG.  22  -  COMPARISON  OF  THE  PRESSURE  DROP  CORRELATIONS 
WITH  THE  DATA  AT  0.65  LB/HR  AIR  RATE 
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CORRELATIONS 


FIG.  25  -  COMPARISON  OF  THE  PRESSURE  DROP 
CORRELATIONS  WITH  THE  DATA  AT  5.10  LB/HR  AIR  RATE 
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FIG.  27  -  COMPARISON  OF  THE  PRESSURE  DROP  CORRELATIONS _ 
WITH  THE  DATA  AT  11.17  LB/HR  AIR  RATE 
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In  comparing  the  predictions  of  Lockhart-Martinel 1 i  and  of 
Kordyban  with  the  data  it  is  seen  that  good  agreement  has  only  been 
obtained  for  the  case  of  0.65  Ib/hr  air  rate.  As  the  air  rate  is  in¬ 
creased,  the  deviations  increase  markedly.  The  data  range  from  30  to 
50  percent  higher  than  predicted  by  Lockhart-Martinel 1 i ,  and  even 
higher  than  by  Kordyban.  Similar  observations  have  been  made  by 
Brigham  (7),  who  reported  pressure  drops  five  to  six  times  those  pre¬ 
dicted  with  flow  inclined  at  12.4  degrees.  Kordyban  (27)  also  found 
that  his  predictions  were  some  20  to  30  percent  higher  than  his  data. 

The  slug  model,  on  the  other  hand,  shows  much  better  agree¬ 
ment  with  the  data.  Although  the  predictions  increase  too  rapidly 
at  the  lower  air  rates,  i.e.,  up  to  40  percent  higher  than  the  data, 
they  return  to  within  25  to  30  percent  above  the  data  over  most  of 
the  range.  The  fact  that  the  predictions  are  consistently  high  makes 
this  model  superior  to  the  correlations  of  Lockhart-Martinel 1 i  and 
Kordyban  for  slug  flow  since  it  provides  an  upper  bound. 

It  must  be  admitted,  however,  that  the  slug  model  only 
approximates  the  actual  slug  flow.  Such  terms  as  slug  acceleration, 
would  tend  to  become  important  with  the  pressure  gradients  that  were 
encountered  in  this  work.  Even  pressure  drop  in  the  gas  phase, 
momentum  loss  in  the  liquid  slugs,  and  others  of  the  simplifying 
assumptions  made  in  the  development  of  the  film  slowdown  problem  could 
become  increasingly  important. 

It  is  unfortunate  that  the  model  could  not  be  derived  from 
purely  mechanistic  considerations.  There  is,  at  present,  no  such 
method  available  for  predicting  film  thickness,  slug  frequency  and  the 
number  of  slugs  in  the  test  section.  As  a  result  it  was  necessary  to 
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employ  actual  data,  or  to  make  approximations  using  the  data,  to  be 
able  to  solve  the  problem.  Nevertheless,  the  model's  predictions  are 
far  more  reliable  and  useful,  under  the  conditions  studied,  than  the 
predictions  of  Lockhart-Marti nel 1 i  and  Kordyban. 


CHAPTER  VII 


CONCLUSIONS 

1.  A  region  of  intermittent  slug  flow,  which  exists  between 
stratified  or  wavy  flow  and  slug  flow,  is  strongly  influenced  by  tube 
i ncl  i nation. 

2.  In  the  transition  from  wavy  or  stratified  flow  to  inter¬ 
mittent  slug  flow  it  is  possible  to  induce  slug  formation  by  artifically 
creating  pressure  fluctuations. 

3.  The  region,  on  a  Baker  plot,  of  regular  slug  flow  is  not 
influenced  by  a  tube  inclination  of  ±7°.  This  indicates  a  stable  flow 
pattern  for  slug  flow. 

4.  Pressure  gradient,  pressure  fluctuations,  slug  velocity,  and 
slug  frequency  are  more  strongly  influenced  by  air  rate  than  by  the 
water  rate,  the  other  remaining  constant. 

5.  Tube  inclination,  within  ±7°,  does  not  have  a  significant 
influence  on  pressure  gradient  (neglecting  static  head),  pressure 
fluctuations,  slug  velocity,  and  slug  frequency.  This  also  shows 
slug  flow  to  be  a  stable  flow  pattern. 

6.  A  sharp  reduction  in  the  slug  frequency  occurs  at  the  transi¬ 
tion  between  plug  and  slug  flow. 

7.  Slug  velocity  and  pressure  gradient  are  influenced  by  the 
transition  from  plug  flow,  through  slug  flow,  to  semi-annular  flow. 

8.  Pressure  fluctuations  in  slug  flow  attain  a  maximum  before 
declining  at  the  onset  of  semi-annular  flow. 

9.  The  observed  slug  velocity  is  up  to  50  percent  higher  than 
the  calculated  total  fluid  velocity. 
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10.  The  Lockhart-Martinel 1 i  and  Kordyban  correlations  predict 
pressure  drop  30  to  50  percent  lower  than  the  slug  flow  data. 

11.  The  pressure  drop  predictions  from  the  proposed  slug  flow 
model  are  25  to  30  percent  higher  than  the  slug  flow  data.  This  pro¬ 
vides  an  upper  bound  for  slug  flow  pressure  drop  predictions. 
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NOMENCLATURE 

A  tube  cross  sectional  area 

D  tube  diameter 

f  Fanning  friction  factor 

F  force 

g  gravitational  acceleration,  gravitational  conversion  factor 

6  mass  velocity  of  gas 

k  thermal  conductivity 

l  film  thickness 

L  length,  mass  velocity  of  liquid 

N  number 

P  pressure 

Q  volumetric  flow  rate 

R  holdup 

Re  Reynolds  number 

r  radius 

S  film  cross  sectional  area 

T  temperature 

t  time 

u  velocity 

w  mass  flow  rate 

X  t(fr)L/(fr)G]1/2  ^ockhart-Martinel 1 1  parameter 

Greek  Letters 
A  difference 

e  liquid  holdup  in  film 
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X  [(pg/0„075)(p^/62.3)]^2  =  parameter  of  Baker  area  plot 

for  two-phase  flow 
y  viscosity 

v  kinematic  viscosity  =  y/p 

tt  3.14159  ... 

p  density 

a  surface  tension 

t  shear  stress 

$  [(fr)TP/(tr)L]1/2  Locl<hart-Martinel  1  i  parameter 

^  (73/a)  [yL(62.2/pL)2]^3  =  parameter  of  Baker  area  plot 

for  two-phase  flow 


Subscripts 
c  constant 

f  from  friction  considerations 

G  gas 

L  liquid 

m  from  momentum  considerations 

s  slug  flow 

TP  two-phase 

t  test 

TOT  total 

1  initial 


2 


fi  nal 
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APPENDIX  -  A 
(1)  AUXILIARY  DATA 


(2)  EXPERIMENTAL  DATA 
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TABLE  A-l 
AUXILIARY  DATA 


PIPELINE  MATERIAL  - 

Tube  Diameter 
Test  Section  Length 
Entrance  Length 
Exit  Length 
Angles  of  Inclination 


CLEAR  EXTRUDED  LUCITE 

=  0.0417  feet 
=  18.29  feet 
=  17.9  feet 
=  4.0  feet 
=  +7°0'  ,  0°0'  ,  -7°0 1 


Viscosity 

Density 


PHYSICAL  PROPERTIES  OF  FLUIDS 

(At  80° F,  1  ATM) 

WATER  AIR  (SATURATED 

WITH  WATER) 

0.900  CP  0.018  CP 

62.20  lb/ft3  0.0708  lb/ft3 


Surface  Tension  67.2  Dynes/Cm 
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TABLE  A-2 
EXPERIMENTAL  DATA 
HORIZONTAL  FLOW 


RATES 

SLUG 

PRESSURES 

AIR 

WATER 

FREQ 

(MIN-1) 

TIME 

VEL. 

AP 

FLUCT. 

(CFM) 

(USGPM) 

(SEC.) 

(FT/SEC) 

(psi) 

(PSI) 

.153 

.53 

72 

5.5 

3.33 

.26 

.12 

.72 

110 

5.0 

3.65 

.38 

.08 

.88 

145 

4.5 

4.06 

.46 

.08 

1.03 

171 

4.3 

4.25 

.58 

.08 

1.17 

201 

4.5 

4.06 

.68 

.06 

1.31 

242 

4.2 

4.35 

.82 

.06 

1.45 

276 

3.7 

4.95 

.94 

.08 

1.58 

305 

3.7 

4.95 

1.06 

.08 

1.70 

335 

3.2 

5.70 

1.18 

.08 

1.81 

361 

3.5 

5.23 

1.30 

.08 

.460 

.53 

41 

2.5 

7.31 

.50 

.40 

.72 

65 

2.2 

8.32 

.72 

.50 

.88 

84 

2.2 

8.32 

.96 

.44 

1.03 

90 

2.0 

9.15 

1.12 

.46 

1.17 

114 

2.0 

9.15 

1.32 

.50 

1.31 

140 

2.0 

9.15 

1.60 

.50 

1.45 

160 

1.9 

9.64 

1.70 

.50 

1.58 

180 

1.8 

10.2 

2.00 

.44 

1.70 

192 

1.7 

10.8 

2.20 

.44 

1 .81 

205 

1.7 

10.8 

2.40 

.50 

.760 

.53 

36 

1.7 

10.8 

.56 

1.20 

.72 

57 

1.6 

11.4 

1.00 

1.20 

.88 

74 

1 .5 

12.2 

1.40 

1.20 

1.03 

85 

1.4 

13.1 

1.80 

1.30 

1.17 

108 

1.3 

14.0 

2.10 

1.20 

1.31 

120 

1.2 

15.2 

2.30 

1.10 

1.45 

144 

1.1 

16.6 

2.60 

1.10 

1.58 

168 

1.1 

16.6 

2.90 

1.00 

1.70 

180 

1.1 

16.6 

3.10 

.90 

1.81 

205 

1.1 

16.6 

3.20 

1.00 

. 
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EXPERIMENTAL  DATA 
HORIZONTAL  FLOW 


RATES 

SLUG 

PRESSURES 

AIR 

WATER 

FREQ 

TIME 

VEL. 

AP 

PLUCT, 

(CFM) 

(USGPM) 

(MIN*1) 

(SEC.) 

(FT/SEC) 

(PSI) 

(PSI) 

1.20 

.53 

16 

.95 

19.3 

.80 

1.20 

.72 

52 

.90 

20.4 

1.50 

1.60 

.88 

60 

.88 

20.8 

1.90 

2.00 

1.03 

79 

.84 

21.8 

2.50 

2.10 

1.17 

102 

.84 

21.8 

3.00 

2.40 

1.31 

115 

.80 

22.9 

3.40 

2.50 

1.45 

132 

.76 

24.1 

3.60 

2.60 

1.58 

152 

.76 

24.1 

4.20 

2.60 

1.70 

162 

.74 

24.7 

4.40 

2.20 

1.81 

206 

.74 

24.7 

4.70 

2.00 

1.82 

.53 

*•. 

.72 

18 

.60 

30.5 

2.00 

2.00 

.88 

50 

.60 

30.5 

2.80 

2.20 

1.03 

77 

.55 

33.3 

3.40 

2.50 

1.17 

100 

.50 

36,6 

4.10 

2.70 

1.31 

115 

.50 

36.6 

4,50 

2.90 

1.45 

132 

.47 

39.0 

5.10 

3.00 

1.58 

150 

.47 

39.0 

5.60 

3.20 

1.70 

160 

.46 

39.7 

5.90 

3.20 

1.81 

174 

.44 

41.5 

6.00 

3.00 

2.63 

.53 

_ 

a. 

■a 

_ 

.72 

- 

- 

- 

- 

- 

.88 

- 

- 

- 

- 

- 

1.03 

52 

.43 

42.5 

4.20 

2.40 

1.17 

80 

.38 

48.1 

5.00 

2.60 

1.31 

108 

.35 

52.2 

5.80 

2.80 

1.45 

120 

.35 

52.2 

6.40 

3.00 

1.58 

- 

- 

- 

- 

- 

1.70 

- 

- 

- 

- 

- 

1.81 

- 

- 

- 

- 

- 
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EXPERIMENTAL  DATA 
POSITIVELY  INCLINED  FLOW  (+7°) 


RATES 

SLUG 

PRESSURES 

AIR 

WATER 

FREQ 

TIME 

VEL. 

AP 

FLUCT. 

(CFM) 

(USGPM) 

(MIN'1) 

(SEC.) 

(FT/SEC) 

(PSI) 

(PSI) 

0.153 

.53 

84 

6.0 

3.05 

.60 

.19 

.72 

116 

4.8 

3.81 

.78 

.14 

.88 

132 

4.4 

4.15 

.82 

.14 

1.03 

180 

4.1 

4.46 

1.06 

.12 

1.17 

220 

4.1 

4.46 

1.20 

.10 

1.31 

240 

3.8 

4.81 

1.34 

.10 

1.45 

276 

3.5 

5.22 

1.50 

.09 

1.58 

302 

3.5 

5.22 

1.70 

.10 

1.70 

327 

3.3 

5.55 

1.80 

.10 

1.81 

346 

3.3 

5.55 

1.90 

.10 

.460 

.53 

56 

2.4 

7.63 

.80 

.54 

.72 

74 

2.3 

7.95 

1.00 

.52 

.88 

95 

1.9 

9.64 

1.24 

.52 

1.03 

134 

1.8 

10.2 

1.50 

.50 

1.17 

168 

1.9 

9.64 

1.70 

.50 

1.31 

188 

1.8 

10.2 

2.00 

.50 

1.45 

212 

1.8 

10.2 

2.20 

.50 

1.58 

240 

1.7 

10.7 

2.30 

.50 

1.70 

260 

1 .65 

11.1 

2.60 

.50 

1.81 

270 

1.6 

11.1 

2.70 

.50 

.760 

.53 

45 

1.5 

12.2 

1.10 

1.10 

.72 

65 

1.4 

13.1 

1.40 

1.20 

.88 

84 

1.35 

13.5 

1.70 

1.20 

1.03 

115 

1.30 

14.0 

2.10 

1.20 

1.17 

132 

1.20 

15.2 

2.40 

1.20 

1.31 

144 

1.15 

15.9 

2.60 

1.20 

1 .45 

160 

1.05 

17.4 

2.80 

1.20 

1.58 

185 

1.05 

17.4 

3.20 

1.20 

1.70 

190 

1.05 

17.4 

3.30 

1.10 

1.81 

200 

1.05 

17.4 

3.60 

1.10 
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EXPERIMENTAL  DATA 
POSITIVELY  INCLINED  FLOW  (+7°) 


RATES 

SLUG 

PRESSURES 

AIR 

WATER 

FREQ 

TIME 

VEL. 

AP 

FLUCT. 

(CFM) 

(USGPM) 

(MIN-1) 

(SEC.) 

(FT/SEC) 

(psi) 

(PSI) 

1.20 

.53 

30 

1.0 

18.3 

1.40 

1.50 

.72 

55 

.90 

20.4 

1.80 

1.80 

.88 

68 

.80 

22.9 

2.30 

2.00 

1.03 

92 

.78 

23.5 

2.80 

2.22 

1.17 

107 

.75 

24.4 

3.20 

2.40 

1 .31 

122 

.75 

24.4 

3.90 

2.50 

1.45 

137 

.75 

24.4 

4.10 

2.50 

1 .58 

153 

.73 

25.1 

4.20 

2.20 

1.70 

166 

.75 

25.1 

4.60 

2.00 

1.81 

180 

.74 

24.7 

4.90 

1.80 

1.82 

.53 

_ 

_ 

_ 

.72 

36 

.62 

29.6 

2.30 

1.80 

.88 

61 

.61 

30.0 

2.80 

2.40 

1.03 

80 

.58 

31.6 

3.50 

2.60 

1.17 

110 

.55 

33.3 

4.30 

2.70 

1.31 

118 

.50 

36.6 

4.80 

3.00 

1.45 

140 

.50 

36 . 6 

5.20 

3.00 

1.58 

151 

.50 

36.6 

5.70 

2.80 

1.70 

168 

.48 

37.4 

- 

- 

1.81 

171 

.47 

36.6 

- 

- 

2.63 

.53 

_ 

_ 

_ 

.72 

- 

- 

- 

- 

- 

.88 

28 

.44 

41.6 

3.80 

2.20 

1.03 

62 

.45 

40.7 

4.60 

2.60 

1.17 

84 

.41 

44.6 

5.20 

2.80 

1.31 

109 

.38 

48.2 

5.90 

2.80 

1.45 

144 

.36 

50.8 

6.50 

2.70 

1.58 

- 

- 

- 

- 

- 

1.70 

- 

- 

- 

- 

- 

1 .81 

- 

- 

- 

- 

- 
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EXPERIMENTAL  DATA 
NEGATIVELY  INCLINED  FLOW  (-7°) 


AIR 

RATES 

WATER 

FREQ 

(MIN’1) 

SLUG 

TIME 

VEL. 

PRESSURES 

AP  FLUCT 

(CFM) 

(USGPM) 

(SEC.) 

(FT/SEC) 

(psi) 

(PSD 

.153 

.53 

* 

.72 

- 

- 

- 

- 

- 

.88 

126 

- 

- 

.06 

.02 

1.03 

160 

- 

- 

.15 

.03 

1.17 

214 

- 

- 

.24 

.03 

1.31 

264 

- 

- 

.33 

.04 

1.45 

316 

- 

- 

.43 

.04 

1.58 

350 

- 

- 

.57 

.05 

1.70 

372 

- 

- 

.65 

.05 

1 .81 

376 

- 

- 

.76 

.06 

.460 

.53 

_ 

_ 

_ 

.72 

50 

2.20 

8.30 

.46 

.28 

.88 

66 

2.00 

9.15 

.60 

.30 

1.03 

115 

1.92 

9.52 

.86 

.32 

1.17 

131 

1.85 

9.90 

1.02 

.38 

1.31 

161 

1.80 

10.15 

1.18 

.40 

1.45 

181 

1.75 

10.4 

1.40 

.40 

1.58 

214 

1.73 

10.6 

1.60 

.40 

1.70 

218 

1.65 

11.1 

1.80 

.40 

1.81 

232 

1.63 

11.2 

2.00 

.40 

.760 

.53 

10 

1.3 

14.1 

<» 

.72 

35 

1.2 

15.2 

.80 

.80 

.88 

58 

1.16 

15.8 

1.30 

1.00 

1.03 

69 

1.10 

16.6 

1.50 

1.20 

1.17 

95 

1.10 

16.6 

1.90 

1.20 

1.31 

105 

1.14 

16.0 

2.20 

1.10 

1.45 

127 

1.15 

15.9 

2.40 

1.00 

1.58 

141 

1.10 

16.6 

2.60 

.90 

1.70 

163 

1.10 

16.6 

2.80 

.90 

1.81 

175 

1.07 

17.1 

3.10 

.80 

B-l 


APPENDIX  -  B 
COMPUTER  PROGRAM 


(1) 

(2) 


CALCULATED  DATA 


D 


2 


C 

C  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  MAKE  PREDICTIONS 

C  OF  PRESSURE  DROP  FOR  HORIZONTAL  SLUG  FLOW  BY  VARIOUS 

C  PRESSURE  CROP  CORRELATIONS 

£ 

C  THE  CORRELATIONS  ARE, 

C  1.  LOCKHART-MARTINELL I 

C  2.  K0R0Y3AN 

C  3.  SLUG  FLOW  MODEL 

C 

£ MAINE  INF  PROGRAM  AND 

C  LQCKH ART-MART  I NELLI  CORRELATION  FOR  TWO-PHASE  FLOW 

C 

COMMON  V  I S C  G , V  I S  C  L , D I  A , P I , G  E  E , D  E  N S  L , D  E  N  S  G , R  E  L  £ , H E  G  P , 

1  KORFUN,  DPL  ,f)PKGR  ,  CO  EFFL  tCOEFFG,  PH  IL ,  PHIG,  RL  ,  N,  M, 

2  RATL ( 25) , RATGC25 ) , I , J , DP SLUG! 60 ) , FREQH { 60 ) ,EREQP (60) , 

_ 3  FREQN(6Q) ,VELH(60)  ,VELP(60  )  ,VELN<60)  , PRES SHI  60)  , _ 

4  PRESSP(6A) ,PRESSN(60) ,K,KK 
LOGICAL  TL,TG,VL,VG 
REAL  M ,  N » K OP  FUN 
C 

C  PHYSICAL  PROPERTIES  AND  COUNTERS  ARE  READ  IN 

C  VISCOSITIES  IN  CENTIPGISE _ 

C  DENSITIES  IN  POUNDS  PER  CUBIC  FOOT 

C  TUBE  DIAMETER  IN  FEET 

RFAD(5,2AA)  V I SCG  »  V [ SCL , DENSG, DE  NSL  *  D 1 A 

R  EAD{  5, 203 )  KK, MM , NN 

C _ 

C  SLUG  FREQUENCIES,  VELOCITIES,  AND  PRESSURE  DROPS  FOR 

C  HORIZONTAL  AND  POSITIVE  AND  NEGATIVE  SLOPES  ARE  READ  IN 

C  THE  LAST  LETTER  BEFORE  THE  SUBSCRIPT  REFERS  TO  THE 

C  TUBE  INCLINATION 

C  H  HORIZONTAL 

C _ P  POSITIVE _ 

C  N  NEGATIVE 

C 

READ { 5. 2041 C i FREQH (K) «V£LH(  K )  , PR E SS H ( K )  )  , K= I , KK ) 

READ! 5, 20  4)  ( {  FREQP(K) ,VELP{K) ,  PRESSPU)  )  ,K=1 ,KK) 

READ ( 5 , 204 ) { (FREQN(K) , V E LN ( K ) , PRE SS N < K ) ) , K= 1 , KK ) 

C _ _ _ 

C  RATES  ARE  READ  IN 

C  WATER  IN  U.S.  G PM 

C  AIR  IN  CJFM 

C 

READ  (5,201)  ( R  ATL ( J ) , J  =  1 »  MM ) 

READ  (5,202  )  (  RA  TG  (  I  )  ,  1  =  1 ,  JMNJ 
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P  1  =  3. 14159 

GEE  =  32 . 17 

c 

c 

c 

HEADINGS  FUR  THE  CALCULATED  DATA  ARE  PRINTED 

WRITE ( 6 ,1000 ) 

WRITE(6,105) 

WRITE  (.6, 1020 

r 

WRITE (6, 104) 

r. 

_ RATES  CONVERTED  TO  1 RS  PER  HOUR 

c 

DO  3  J=  1 ,  MM 

3  PATL ( J ) =  R  A  T  L ( J )*6^. 0*.  1336*DENSL 

DO  2  1=  1 , NN 

2  RATGI I ) =RATG( I ) *60.Q*DENSG 

K=1   

DO  1  1=1, NN 

DO  1  J  = 1 , V  M 

TL=. FALSE. 

TG=. FALSE. 

VL=. FALSE  . 

VG=. FALSE. 

c 

c 

c 

THE  REYNOLDS  NUMBERS  FOR  THE  GAS  AND  LIQUID 
ARE  CALCULATED,  ASSUMING  EACH  PHASE  TO  FLOW 

PHASES 

ALONE 

c 

REGP={4.0*RATGI  I)  ) / < VI SCG*D I A*P I *2 .4191 ) 
RFLP= (4.0*RATL ( J)  )  /  (  V  ISCL*D  IA*P 1*2.4191) 

c 

c 

c 

DETERMINATION  WHETHER  EACH  CF  THE  PHASES  IS 
TURBULENT  OR  LAM  I NAR I  V I  SCO US)  FLOW 

IN  EITHER 

c 

I F ( R  E  GP . L  E . 2  000 . 0 )  VG=.TRUE. 

IFtRELP.LE. 2000.0)  VL=.TRUE. 

c 

IF (RELP.GT. 2000.0)  TL  = . TRUE . 

IFtREGP.GT. 2000.0)  TG=.TRUE. 

c 

c 

c 

DEPENDING  ON  THE  CHOICE  OF  FLOW  MECHAN I SM, ONE  OF  THE 

FOUR  SUBROUTINES  IS  CALLED  TO  CALCULATE  THE  PARAMETERS 

IF(TL.AND.TG)  CALL  TURTUR 

IF (TL  .  AND. VG)  CALL  TURVIS 

IF( VL.AND.TG)  CALL  VISTUR 

r 

IF ( VL . AND.VG)  CALL  VISVIS 

L 

c 

FRICTION  FACTORS  ARE  CALCULATED  BY  THE  BLASIUS  EQUATION 

'  .  = 

.  i 

r .  . 
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C 

FFG=COEFFG/ ( REGP**M ) 

FFL=COEFFL/(  RELP**N » 

C 

£ _ SINGLE  PHASF  PRESSURE  ORGP  CALCULATIONS _ 

C  BY  FANNING  EQUATION 

C 

QPL=  (  32 .0*FFL*(  RATL  (  J  )**2i  )  /  (  DEJiSL*  (  D  I  A**  5 )  *(  PI**2) 

I  *GEF*(  3600 ,0**2)  *144.0) 

DPG=(32 .0*FFG*( RATG ( I )**2) ) / ( DENSG* ( D I A**5 ) * { P I **2 ) 
_ 1  *GEF*(3600.0**2)*144.C ) _ 

C 

C  TWO-PHASE  PRESSURE  DROP  CALCULATIONS  FOR  LOCKHART- 

C  MART! NELL  I  CORRELATION 

C 

DP  TP  L- < PH IL**2 ) *D  PL 

_ DPTPG= ( PH IG**2 ) *UPG _ 

C 

C  THE  KCRDY BAN  CORRELATION  IS  CALLED 

C 

CALL  KCRDY 
C 

£ _ THF  FXPFR  TMFNTAL  DATA  FOR  OVERALL  PRESSURE  DROP 

C  ARE  CHANGED  TO  PRESSURE  GRADIENT  PSI/FT 

C 

PRESSHI K) =PRESSH(K )  /  i  8 . 2  9 
PRESSP( Kl -PRESS P( K ) / 1 8*  29 
PRES  SN ( K )=PRESSN{K)/18.29 

C _ 

C  THE  SLUG  MODEL  IS  CALLED 

C 

CALL  SLUG 

IF(( I.EQ.4) .ANO.t J.EQ.l))  GO  TO  4 
GO  TO  5 

4  WRITE  (6, 1"0^) _ _ _ 

WRITE (6, 104) 

5  CONTINUE 

WRIT  E (6  > 1 03)  RATH  Jit RATG (I ) ,PPTPL,DPKOR,PPSLUG(K) t 
1  PRESSHCK) tPRESSP(K) *PRESSN(K) 

K  =  K  +  1 

_ LK=KK+ 1 _ _  _  _  _ 

IF(K.EQ.LK)  GO  TO  10 
1  CONTINUE 

102  FORMAT (34X» "CALCULATED  DATA1//) 

103  FORMAT  (2XfF8.3~,  F10.3»6F10*5) 

104  FORMAT (5X, 'RATES  ( LB/HR  )*, 1 5X, ' PRESSURE  GRADIENT  (PSI*, 

l »  / FT)  *  /5Xt  »  WATER  AIR*  ,  6X,  «  L-M_ KORDY_ S L U  •  ± 
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2 ' G  HORTZ  POS  NEG1/) 

105  FORMAT (  3  7  X  f  ’  T ABLE  8-1'/) 

200  FORMAT { 5F10. 5 ) 

201  FORMAT (4F15.5/4F15.5/2F15.5 ) 

202  FORMAT (4F15«5/2F15»5) _ 

203  FORMAT (313) 

20 A  FORMAT ( 3 F 1 0 *  2 ) 
lOOO  FORMAT ( 1H 1 ) 

10  STOP 
END 
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SUBROUTINE  VISVIS 
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C 

C  SUBROUTINE  VI SCOUS- VISCOUS 

C  CALCULATION  OF  THE  PARAMETERS  OF  THE  LOCKHART -MAR T I N- 

C _ ELLI  CORRELATION  FOR  LIQUID  AND  GAS  BOTH  IN  VISCOUS 

C  FLOW 

C 

COMMON  VI SCG, VISCLfOI A, PI , G E E , DENSL , OENSG , R EL P , REGP , 

1  KORFUN f  DPL  > DPKQR  » CO EF FL , COE FF G , PH  I L , P H 1 G , RL , N , H , 

2  P  ATL ( 2  5 )  ,RATG(25),  I , J , DPSLUG ( 6H ) tFREQH(60)  »FREQP(60)  , 

_ 3  FREQNI 60)  tVELH(6C  )  »  V  E  L  P  (  6  0  )  , V ELN ( 60 )  ,  P R F S SH ( 6^ )  , _ 

4  PRESSPC 60) ,PRESSNI60) ,K,KK 
REAL  M » N 
C 

C  8LASIUS  EQUATION  NUMERATORS  AND  EXPONENTS 

C 

_ C0EFFL  =  l6.r> _ 

COEFFG=16.0 
N  =  1 .0 
M- 1 »  G 
C 

C  CALCULATION  OF  PARAMETERS  X  »  RL  »  AND  PHIL  OF  THE 

C_ LOCKHART -MART  I NELL  I  CORRELATION 

C 

X-SQRT ( ( RATL { J ) /RATGU )  )*{ D ENSG/DENSL )*IVISCL/VISCG)) 

X  =  A  L  0  G 1  0  (  X  ) 

PHI L= .42  5  3-. 42 2  7*X+. 21 18*1 X**2)-. 3480E-01*( X**3 ) 

1  -.6754E-C2*(X**4) 

RL  =  -  .  64  31  +  . 4816*X-  .  12  89  *{  X**2  )  +_.  2323E-0  1»(  X**3  ) 

1  2267E-n2* { X  *  44 ) 

C 

C  ANTILOGS  OF  XT  RL t  AND  PHIL 

C 

X=10 .o**x 

_ RL=10. 9**RL _  _ 

PHIL=10.9**PHIL 

PHIG=X*PHIL 

RETURN 

END 
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SUBROUTINE  VISTUR 
C 

C  SUBROUTINE  VISCOUS-TURBULENT 

C  CALCULATION  OF  THE  PARAMETERS  OF  THE  LOC KHA RT -MAR T I N- 

C _ ELL  I  CORRELATION  FOR  THE  LIQUID  IN  VISCOUS  FLOW  AND 

C  THE  GAS  IN  TURBULENT  FLOW 

C 

COMMON  VI SCG, VISCL ,  D  1  A , P  I  , G E E , DENSE , OEN SG , R ELP , R EGP , 

1.  KQRFU  N  *  D  PL , DPK  OR  f COEFFL , CGEFFG,  PHIL,  PHIu,RL,N,M, 

2  RATL ( 25 ) ,RATG(25) , I , J , DPSLUG < 60 ) , F REQH ( 60 ) , FREQP ( 60 ) , 

_ 3  FREQN  (60  )  «  V  ELH { 60  )  ,VELP  (60  )  ,VELN(6Q)  .PRESSHtoO)  , _ 

4  PRESSPI60) ,PRESSN(6Q) ,K,KK 
REAL  M , N 
C  _ _ 

C  BLASIUS  EQUATION  NUMERATORS  AND  EXPONENTS 

C 

_ CQEFFL= 16.0 _ 

C0EFFG=0 .046 
N  =  1 .0 
M=0. 2 
C 

C  CALCULATION  OF  PARAMETERS  X , RL  »  AND  PHIL  OF  THE 

C  LOCKHART-MART INELLI  CORRELATION _ 

C 

X=SQRT I { COEFFL/COEFFG) * (REGP**< -0.8) ) *{ RATL  < J I/RATGi I ) ) 
1  *( DENSG/DENSU *< VI SCL/VISCG ) ) 

X  =  A  L  OG 1 0  I  X ) 

P  H I  L  =  .  5396-.4661*X+.  153  L  *(  X**2 )-.  2207E-01*{  X  ^  3  ) 

_ 1  +. 2642F-"3*<  X**4) _ _  _ 

RL=-.6431+. 481 8*X-. 1289*{ X**2) + . 2323E-0 1* ( X**3) 

1  -.2267E-C2*{ X*^4 ) 

C 

C  ANT  I  LOGS  OF  X,  RL,  AND  PHIL 

C 

_ X=  10 ,C**X  _  _ 

RL=10.0**RL 
PH  I L= 10 • 0**PH I L 
PHIG=X*PH IL 
RETURN 
END 
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C 

SUBROUTINE  TURVIS 

c 

c 

r 

SUBROUTINE  TURBULENT-VISCOUS 

CALCULATION  OF  THE  PARAMETERS  OF  THE  LOCKHART -MART  I N- 
ELL I  CORRELATION  FOR  THE  LIQUID  IN  TURBULENT  FLOW  AND 

c 

r 

THE  GAS  IN  VISCOUS  FLOW 

COMMON  VI  SCG,  VISCL  ,DIA,  P  I  ♦  G  EE  ,  DF.NSL  .  DENSG  »  R  EL  P  *  REGP  , 

1  KORFUN»  DPL  »  DPKGR  , CO EF FL , COEFFG, PH IL , PH I G , RL , N , M , 

2  RATL (25)  ,RATG( 25) ,  I  , J,DPSLUG<60 ) ,FREQH(60) ,FREQP(69)  , 

3  FREQN(6^  )  , VELHC60  )  . VELP<60  )  * VELN(60)  »PRESSH( 60) i 

C 

4  PRESSP(60) ,PRESSN(60) ,K,KK 

REAL  M,N 

C 

c 

BL AS  I  US  EQUATION  NUMERATORS  AND  EXPONENTS 

COEF  F  L  =  Q .346 

CO E F  FG= 1 6  •  0 

N=0 . 2 

M=  1  •  0 

c 

c 

c 

CALCULATION  OF  PARAMETERS  X  t  RL  *  AND  PHIL  OF  THE 

1 OCKHART -MART INELL I  CORRELATION 

c 

X=  SQRT ( ( COEFFL/ COEFEG ) * ( RELP**0 . 8 3 * ( R ATL ( J ) /RATGt  1  I  )* 

I  ( DENSG/DENSU*< VI SCL/V ISCG 3 ) 

X  =  A  LOG 1 0 ( X) 

PHIL=.5433-.4406*X+. 1443* ( X**2)-. 2778 E-G1*(X**3> 

1  +. 1663E-Q2*( X**4 ) 

c 

RL=-. 64  31 +  . 4818* X-. 12  89* ( X* *2) +  .2  32 3E-01*(  X**3) 

1  2267E-02*( X**4) 

c. 

c 

ANTILGGS  OF  X,  RL ,  AND  PHIL 

X=10 ,c**x 

RL= 10.0**RL 

PHI  L  =  1 A  ,Q**PHIL 

PHI  GJ=X*PH  I  L 

RETURN 

END 
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SUBROUTINE  TURTUR 
C 

C  SUBROUTINE  TURBULENT-TURBULENT 

C  CALCULATION  OF  THE  PARAMETERS  OF  THE  LOCKHART -MAR T I N- 

C_ ELLI  CORRELATION  FOR  BOTH  THE  LIQUID  AND  THE  GAS  IN 


C 

r 

TURBULENT  FLOW 

COMMON  VI  $CG»  VISCL  ,DI  A,P  I  ,GE£»  DE.NSL  .DENSG, R.ELP*.REGP» 

1  KORFUNt DPL ,DPKOR  , COEFFL , COEFFG, PH IL , PH  I G* RL *N , 

2  RATL  (25)  ,RATb(25  )  ,  I  ,  J,  DPS  LUG  (60  )  f  FREQH(60)  ,FRE(jP(60)  , 
_ 3  FREQN ( 60 ) , VELH160)  ,VELP(60  )  ,VELN<6C )  ,PRESSH(6G}  , 

C 

4  PRESSPI 60) ,PRESSN{60) ,K,KK 

REAL  M »  N 

C 

r 

BLASIUS  EQUATION  NUMERATORS  AND  EXPONENTS 

COEFFL=0  .  n4  6 

COE  F  FG=  T .  °46 

N=0 . 2 

M=n .  2 

c 

C 

C 

CALCULATION  OF  PARAMETERS  X  » RL  »  AND  PHIL  OF  THE 

LOCKHART -MART  IN  ELI.  I  CORRELATION 

c 

X  =  SQRT ( ( (RATL( J )/RATG ( I  )  )**1.8 ) *{ DENSG/DENSL ) 

1  ( VISCL/V  ISCG  )**0.2 ) ) 

X= ALOGIO (X) 

PH  I L=. 62 7 5-. 5  01 3*X+.  1308*1 X**2 ) - .  8409 E-02*  (  X**3  i 

1  -.2251E-G2*! X**4) 

C 

RL=- . 643 1+ . 4 8 1 8* X- .  1 2  89* ( X* * 2  >  +  . 232 3E -0  I’M  X** 3 ) 

1  -.2267E-Q2JM  X**4) 

C 

c 

ANTILOGS  OF  X,  RL»  AND  PHIL 

X= 10* o**x 

RL=10.0**RL 

PHIL=10.0**PHIL 

PHI G= X*  PH  I  L 

RETURN 

END 
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c 

SUBROUTINE  KURD Y 

c 

c 

c 

KOROYBAN  CORRELATION 

COMMON  VI  SCO, VI SCL, 01  A, PI ,GEE, DENSL , DENSG,RELP ,REGP , 

1  KORFUN  »  OPL  »  DPKCR  , COEFFL , COEFFG, PH I L , PHI G , RL , N, M, 

2  RATL  (25)  ,.RATG  I  25  )  ,  I  ,  J  ,  DPSLUG  I  60  )  »  F  REQH  L60  )  ,  FREQP  (  60.)  , 

3  FREQNI6Q) , VELHI60) , VELP(60) ,VELN(60) , PRES SHI 60) , 

4  PRESSPI60)  ,PRESSN( 60)  ,K,KK 
REAL  M «  N « K  CP  F  UN 


C 

C  CALCULATION  OF  THE  KQRDYBAN  FUNCTION 


C 

C 

KQRFUN= I .0+( DENSL/DENSG  )  *(RATG(  I) /RATL( J)  ) 

THE  SINGLE  PHASE  LIJUIO  PRESSURE  DROP  FROM  THE 

C 

LOCKHART-MARTINELLI  CORRELATION  IS 

USED  TO  CALCULATE 

c 

c 

THE  PREDICTED  PRESSURE  DROP  BY  THE 

KORDYBAN  CORRELATION 

c 

r 

KOROYBAN  TWO  PHASE  PRESSURE  DROP 

DPKOR  =DPL*<  KORFUN**.  75).  ... 

RETURN 

END 
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C 

C 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 


c 

c 


c 

c 

c 

c 

c 


SUBROUTINE  SLUG 

CALCULATION  OF  THE  PREDICTED  PRESSURE  GRADIENT  BY  THE 
SLUG  FLOW  MODEL 


COMMON  VI SCGt VI SCL ,DI A, PI , G EES DENSL , DENSG » RELP , REGP , 

1  KO  RFUN  rOPL, DPKOR  , COE F FL , COEFFG , PH  I L , PH I G , RL , N , M , 

2  RATL ( 23)  ,RATG( 2  5 ) ,  I , J , DP SLUG (60 )  , FREQH<6P )  ,F  REQP ( 60)  , 

3  FREQNf 60) , VELH ( 60 ) ,V£LP(60 ) ,VELN(60) ,PRESSH(60> , 

4  PRESSPI 60) , PRESSNI60) ,K»KK 

DIMENSION  RESl  UG( 60) . FFSLUG<  60 )  .DPLSL (60 ) , DPMQMH ( 60 )  , 

1  QL ( 60) »  QG( 60 ) t  V {  60) t  FR ACTL ( 60 ) »  NSLUG ( 60 ) 

REAL  M,N 
REAL  NS  LUG 

IFCFREQHCK) .EQ.0.0 J  GG  TO  2 

THE  MASS  RATES  ARE  CONVERTED  TO  VOLUMETRIC  FLOW  RATES 
CUBIC  FEET  PER  SECOND 

QL  (  J ) =  R A T L ( J ) / ( DE NS L *3 600.0  ) 

QG< I )=RATG( I )/( 0ENSG*36Q0.Q ) 

THE  FRACTION  LIQUID  I N  THE  PIPE  IS  CALCULATED 


FRACTLI K ) =QL { J ) / i QL ( J )+QG( I ) ) 
OVERALL  FLUID  VELOCITY 


V  iK)=4.n*(QLiJi^wG(I  )  )  /  <  P  I  *(DI A**2)  1 


THE  REYNOLDS  NUMBER  IS  CALCULATED  AND  USED  TO  CALCULATE 
THE  FRICTION  FACTOR  BY  THE  BLASIUS  EQUATION 
(TURBULENT  FLOW  IS  ASSUMED) 

RESLUG(K)=DIA*V(K) »DENSL/ i VISCL *6.72E-Q4) _ 

FFSLUG( K) =0.0 791/ (RESLUG (K) **0.2 5) 

PRESSURE  DROP  FROM  LIQUID  AT  FLUID  VELOCITY 


DPLS  L ( K ) =  2 . ^  *FFSL UG ( K ) *D EN  SL * ( V ( K ) **2 1 / <  GEE *D I  A* 1 44 . 0 ) 
FREQH(K)  =  FREQH(K)/60 . 0 _ 

THE  NUMBER  OF  SLUGS  IN  THE  TEST  SECTION  IS  CALCULATED 

NSLUG no  = F p  E QH I K ) / V ( K ) *  1 8 . 29 

CALCULAT I  ON  OF  THE  EXPONENT  I AL  TERM  FOR  THE  FILM _ 
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B  -  12 


C 

C 


L 

C 

C 


C 

C 

c 


c 

c 

c 


SLOWDOWN 


B=<P|**2 ) *{ VI$CL*6.72*1QE-47DENSL)/(FREQH(K) *4.0*( C  (RL 
l-FRACTLC  K)  )  *DI  A  )**2  )  ) 


SUMMATION  IN  MOMENTUM  TERM 

SUM  1=^.0 
DO  1  1 1  =  1 f  5 

A= ( 2.0*FL0AT< I  I  )-l  .0)**2 
SUM1  =  SIJM1  +  1  .0/A*FXP(-A*B) 


1  CONTINUE 


CALCULATION  OE  THE  MOMENTUM  LOSS  OE  THE  FILM 


DPMOMH ( K ) =NSLUG ( K )  *DENSL* ( RL-FRACTL { K ) ) *( V( K )**2 ) / 
1  <GEE*144.0*18.?9]*<l.?-64.0/(P 1**4 ) * [ SUM1**2 JJ _ 


TOTAL  SLUG  PRESSURE  DROP  PREDICTION 


DPSLUG ( K  )  =  OPLSL(K)*FKACTL (K)+OPMGMH(  K) 
GO  TO  3 

2  DPSLUG (K) =0.0 


3  RETURN 
END 
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APPENDIX  -  C 
SAMPLE  CALCULATIONS 

(1)  LOCKHART-MARTINELLI  CORRELATION 

(2)  KORDYBAN  CORRELATION 

(3) 


SLUG  FLOW  MODEL 


I 
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PROBLEM:  To  predict  the  two-phase  pressure  drop  in  a  0.5  inch  diameter 
tube. 

Water  Rate  =1.31  USGPM 
Air  Rate  =0.76  CFM 

PHYSICAL  PROPERTIES 
(AT  80° F  and  1  ATM) 


VISCOSITY 

DENSITY 

(CENTI POISE ) 

(LB/FT3) 

Air 

0.018 

0.0708 

Water 

0.900 

62.20 

1.  Solution  by  Lockhart-Martinel 1 i  Correlation 
Convert  rates  to  Ib/hr 
WL  =  Ql  (60.0) ( . 1 336)  pL 

=  (1 .31 ) (60.0) ( . 1336) (62.2) 

=  653.16  Ib/hr 
WG  =  Qg  (60.0)  pG 

=  (0.76) (60.0) (0.0708) 

=  3.228  Ib/hr 

Calculation  of  the  gas  and  liquid  Reynold's  Numbers  if  each  is  assumed 
to  flow  alone. 

ReL  =  (4.0)  Wl/p(2.4191)ttD 

=  (4.0)  (653. 16)/ (0.900)  (2.4191  )tt(0.0417) 

=  9160.0 

ReQ  =  (4.0)  WQ/y (2.4191  )ttD 

=  (4.0)  (3.228)/ (0.018)  (2.4191  )tt(  .041  7 ) 

=  2263.8 

Since  the  Reynolds  numbers  are  above  2000,  both  phases  are  in  turbulent 


flow. 


The  Blasius  numerators  and  exponents  for  this  case  are: 


liquid  gas 

n  =  0.2  m  =  1 .0 
CL  =  0.046  CG  =  0.046 
fL  =  0.046/(ReL)2’0 
=  0. 046/(9160. 0)2 
=  0.0074 


Pressure  drop  in  the  liquid  phase 

APl  =  (32.0)  fLWL2/PLD57T2gc  (3600.0)2(144.0) 

=  (32. 0)  (. 0074)  (653. 16)2/ (62. 2)  (0.041 7)  V(  32. 17) 

(3600.0)2(144.0) 

=  0.02180  psi/ft 


Calculation  of  Parameter  X  = 


( 


(— ) 
^al;l 

7MET 

^AL;G 


) 


0.5 


for  turbulent-turbulent  flow 


W 


^,^■^^0.2 


,653.16x1.8  ,0.0708x, 0.900x0. 2 
'  ^3 . 228  ’  y  62.2  M0.018; 

X2  =  35.3 


X  =  5.94 

from  Lockhart-Marti nel 1 i ' s  Coordinates  of  $  and  versus  parameter  X 

(curve-fit  in  the  Computer  Program  in  Appendix  B) 

/  AP  x 


. 
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R^  =  liquid  holdup  =  0.4596 


Therefore 


APtp  -  $L  APl 


=  (2.057;  (0.02180) 

=  0.09226  psi/ft 

Pressure  Drop  Prediction  by  Lockhart-Martinel  1  i  =  0.09226  psi/ft. 

2.  Kordyban  Correlation 
Calculation  of  the  Kordyban  function 
K0RFUN  =  1.0  +  (PL/pg)(WG/WL) 

=  1.0  +  (62. 2/0. 0708)(3. 228/653. 16) 

=  5.342 

Calculation  of  single  phase  (liquid)  pressure  drop.  Kordyban  assumed 
only  turbulent  flow  in  smooth  tubes  having  a  Reynold's  number  of  less 
than  100,000. 


Therefore  Re  = 


0.079 


Re 


0.25 


To  extend  the  range  of  Kordyban 's  correlation,  the  Blasius 
equation  used  by  Lockart-Marti nel 1 i  was  assumed  to  be  more  applicable 
Therefore,  AP  from  Lockhart-Martinel 1 i  is 
AP,  =  0.02180  psi/ft. 

Therefore 

,0.75 
,0.75 


AP-j-p  =  APl  (K0RFUN)' 


=  0.02180  (5.342) 

=  0.07659  psi/ft 

Pressure  Drop  Prediction  by  Kordyban's  Correlation 


=  0.07659  psi/ft 


3.  Slug  Model  Theory 

GIVEN:  Slug  frequency  =  1 20/Mi n 

Calculation  of  the  fraction  liquid  in  the  total  flow. 

FRACTL  =  Ql/(Ql  +  Qg) 

/  10.50  \ 

MO.  50  +  45.6' 

=  0.187 

Calculation  of  fluid  velocity 

ttD2 

U  =  (Ql  +  Qq)/  ^ 

=  (56 . 1 0)  (4 . 0)/tt(  .041 7)2(3600) 

=  11.4  ft/sec 

Calculation  of  frictional  pressure  drop  term. 

Fluid  Reynolds  number,  assuming  fluid  density  to  be  that  of  the  liquid. 
Re$  =  DupL/y 

=  (0 . 041 7) (11 .4) (62 . 2)/ (0 .900) (6.72  x  10"4) 

=  48935.1 

Friction  factor  by  Blasius  equation 
f  =  0. 0791/Re0,25 

=  0. 0791/(48935. I)0,25 
=  0.00532 

Frictional  pressure  drop  of  the  slugs 

o  Q, 

AP$  =  ( ( 2 . 0 ) f pL  uVgcD(144.0))(Q  $  -Q  ) 

= ((2. 0) (0. 00532) (62. 2) ( 11. 4)2/ (32. 17) (0.041 7) 

( 1 44 . 0 ))( 0 . 1 87 ) 


=  0.0834  ps i /ft 
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Calculation  of  momentum  loss  term.  Number  of  slugs  in  the  test  section 

N  =  slug  frequency  L 
s  u  Lt 

=  (if^)08.29) 

=  3.20  slugs/section 

Calculation  of  the  exponential  term  for  film  slowdown 

B  =  (tt^(6.72  x  10_4)/Pl)/(  (Slug  frequency)  (4.0) 

(Rl  -  FRACTL)D2) 

=  (tt2(0.900)(6.72  x  10"4)/(62.4)  )/((2.0)  (4.0)  ( .4596  -  .187) 

(0.0417)2) 

=  0.9296 

Calculation  of  summation  in  the  momentum  term 

oo  p 

SUM  =  l  (1.0/(2n-l)2)e"^2n“^  B 
1 

Since  this  term  converges  very  quickly,  the  summation  was  taken  from 
1  to  5,  only 

SUM  =  0.3947 

Therefore , 

Total  momentum  loss  term 

APm  =  (NspL^RL  "  FRACTL)u2/gc(144.0)Lt)((1.0  -  (64. 0/tt4)  )  (SUM)2) 

=  ( (3.20) (62.2) (0.4596  -  0. 1 87) ( 1 1 . 4) 2/ (32 . 1 7) ( 1 44. 0) (1 8. 29) ) 
((1.0  -  (64. 0/tt4)  )  (0. 3947) 2 ) 

=  0.07493  psi/ft 

Therefore  the  total  pressure  drop  by  the  slug  flow  model  is 

APTnT  =  AP  +  AP 
TOT  s  m 

=  0.0834  +  0.0749 
=  0.1583  psi/ft 


■ 


Actual  Pressure  drop  was  0.1258  psi/ft 


Relationshi p 

Lockhart-Martinel 1 i 

Kordyban 

Slug  Model 


Comparison 


0.09226 

0.07659 

0.1583 


Deviation  from 
actual 

-26.6% 

-39.1% 

+25.8% 


Actual 


0.1258 
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APPENDIX  -  D 
INSTRUMENT  CALIBRATION 

(1)  ROTAMETERS 

(2)  PRESSURE  TRANSDUCERS 
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ROTAMETER  CALIBRATION 

Factory  calibration  curves  were  supplied  with  each  of  the 
Brooks  Rotameters 


AIR  WATER 

#6709-31235  #6709-31236 

#6709-2288-1 

No  curves  were  available  for  #624-675-1. 

All  of  the  meters  were  recalibrated  for  the  respective  fluids 

at  80° F. 

Calibration  of  the  air  Rotameters  was  carried  out  with  the 
use  of  a  V0L-0-FL0  area  meter  (National  Instrument  Laboratories ,  Inc.) 
(Mod.  20  ER)  which  has  a  range  of  zero  to  300  scfh.  This  meter  was 
previously  calibrated  with  a  volume  displacement  meter. 

The  results  are  shown  in  Table  D-l  and  are  plotted  in  Figures 

28  and  29. 

Calibration  of  the  water  Rotameters  was  carried  out  by  measuring 
the  flow  into  a  graduated  cylinder  over  specific  time  intervals.  The 
flow  rates,  in  milliliters  per  minute  were  converted  to  U.S.  gallons  per 
minute  versus  the  scale  readings.  The  resultant  data  are  shown  in 
Table  D-2  and  are  plotted  in  Figures  30  and  31. 
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TABLE  D-l 
ROTAMETER 
CALIBRATION  DATA 

WITH  WATER  SATURATED  AIR  AT  80°F 


Brooks  Rotameter  #6709-2288-1 

Model  8-1100,  Tube  No.  R-8M-25-4,  Float  8RV3 


Scale  Reading 


cfm 

(80°F,  14.65  psia) 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 


.392 

.650 

.917 

1.20 

1.475 

1.82 

2.18 

2.635 

3.15 

3.70 


Brooks  Rotameter  #6709-31235 

Model  1110  Tube  R-6L-600-G1 ,  Float  316SS  (RS  Spool) 


Scale  Reading 


cfm 

(80°F ,  14.65  psia) 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 


.1525 

.249 

.317 

.385 

.460 

.530 

.600 

.670 

.765 

.867 


0NI113S  y313W\/10B 


CALIBRATION  OF  BROOKS  ROTAMETER  #6709-31235 


MODEL  NO.  8-1100,  TUBE  NO.  R-8M-25-4 
FLOAT  8RV3 


CALIBRATION  OF  BROOKS  ROTAMETER  #6709-2288-1 
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TABLE  D-2 

ROTAMETER 

CALIBRATION  DATA 
WITH  WATER  AT  80° F 

Brooks  Rotameter  #624-675-1 

Model  8-1100,  Tube  No.  R-8M-25-4,  Float  8-RS-14 


Scale  Reading 

USGPM 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

.526 

.722 

.887 

1.024 

1.16 

1.305 

1.445 

1.58 

1.70 

1.80 

Brooks  Rotameter  #6709-31236 

Model  8-1110,  Tube  No.  R-6L-600-G1 ,  Float  316SS(Rs  Spool) 


Scale  Reading 

USGPM 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

.0395 

.0534 

.077 

.099 

.116 

.133 

.152 

.1695 

.1915 

.212 

MODEL  8-1110  TUBE  NO.  R-6L-600-G1 
FLOAT  316SS  (RS  SPOOL) 

WATER  TEMPERATURE  80° F 
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WATER  FLOWRATE  (U.S.  GALLONS  PER  MINUTE) 

F!G.  31  -  CALIBRATION  OF  BROOKS  ROTAMETER  #624-675-1 


PRESSURE  TRANSDUCER  CALIBRATION 


The  two  Statham  Pressure  Transducers  Model  80TC  ±15  -  350 
#10153  and  #3830  were  calibrated  in  the  following  manner.  Pressure 
was  applied  to  the  positive  side,  with  the  negative  side  open  to  the 
atmosphere,  from  the  compressed  air  system  and  was  measured  with  a 
mercury  U-tube  manometer.  The  output  was  recorded  on  a  Hewlett- 
Packard  (Moseley  Autograf  Model  7100B)  stripchart  recorder.  The  re¬ 
sults  are  shown  in  Table  D-3  and  are  plotted  in  Figure  32. 
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TABLE  D-3 
CALIBRATION  DATA 
FOR 

STATHAM  PRESSURE  TRANSDUCERS 
MODEL  80  TC  ±15  -  350 
#10153  AND  #3830 

Hewlett-Packard  Carrier  Preamp  Setting:  Atten  =  20 
Hewlett-Packard  Autograf  Recorder  Setting:  5 V 


Autograf  Reading 
1.0 
2.0 
3.0 
4.0 
5.0 
6.0 
7.0 
8.0 
9.0 


Pressure  (PSIG) 
1.03 
2.06 
3.10 
4.12 
5.15 

6.19 

7.20 
8.24 
9.27 


H-P  AUTOGRAF  READING 


